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Turbulent Heat and Mass Transfer from 


Stationary Particles’ 


I. S. PASTERNAK? and W. H. GAUVIN® 


Convective heat and mass transfer rates for 20 
shapes suspended in various orientations in a hot 
turbulent air stream were correlated by the equation: 


j'p = 0.692 (Re’’)~°-486 for 500 < Re’’ < 5000 


for turbulent intensities between 9 and 10%, using a 
new characteristic dimension, L”, in the Reynolds 
Number. This dimension, defined as the total surface 
area divided by the maximum perimeter perpendic- 
ular to flow, was qualitatively explained from bound- 
ary layer considerations and was found to correlate 
the data for different shapes available in the liter- 
ature. Qualitative pictures for the relative rates of 
heat and mass transfer around the various shapes 
were also obtained. 


Ale: amount of experimental evidence has been reported 
in the literature on heat and mass transfer coefficients for 
standard bodies such as cylinders, spheres and flat plates, but 
little information is available on other shapes such as disks, 
prisms, rings, granular particles or fibres. Since the latter are 
finding increasing use in practical engineeri ; applications, such 
as flash drying and many solids-gas contacting operations, modi- 
fications to existing correlations are required to account for the 
influence of the shape on the transfer phenomena. A more funda- 
mental approach to the problem would be to define a uniform 
particle characteristic dimension by means of which all data on 
heat and mass transfer could be correlated, irrespective of the 
shape of the particle. 

A survey of the published literature shows that many 
“particle diameters” have been proposed by various authors in 
an attempt to define a particle of non-standard shape“: * * ”), 
The most commonly used diameter for a single particle (and the 
easiest to determine) appears to be the equivalent diameter, 
d,, the diameter of a sphere of equal volume, while d,, the dia- 
meter of a sphere of equivalent surface area, has also been used. 
In the field of multiparticle technology, such as in the flow of 
fluids through packed beds, Wadell’s sphericity y‘, Carman’s 
$; (= W) ‘© as well as Leva’s \ (= 1/) have become popular 
as a means of designating the average shape of a particle. How- 
ever, these particle diameters and shape factors all suffer from 
a common disadvantage —- they constitute a fixed property of 
a particular shape without reference to the particle orientation. 


\Manuscript received October 29, 1959; Accepted January 11, 1960. 
2Imperial Oil Limited, Sarnia, Ont. 

3McGill University and Pulp and Paper Research Institute of Canada, 
Montreal, Que. 
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In 1940, Powell showed that both particle shape and orien- 
tation had a large effect on the rate of heat and mass transfer, 
but he made no attempt to correlate either of these factors into 
the existing transfer equations i.e. the Nusselt equation®: ' ™ 
and the Chilton-Colburn equation “?: 1%. 19, “The validity of 
these equations is now well established and they are commonly 
used to correlate the individual heat and mass transfer data for 
single spheres, cylinders and flat plates, provided that the dia- 
meter, for both the sphere and the cylinder, or the length of the 
flat plate in the flow direction, is used as the characteristic di- 
mension in the Reynolds Number (and in the Nusselt Number 
if it is employed). 

In 1953, Krischer“® suggested the use of a characteristic 
length which would take into account both the shape and orienta- 
tion of a particle and defined a “flow length”, Lx, as the length 
of the path covered by the streaming medium along the transfer 
surface. Using this flow length in the Nusselt equation, he was 
partially successful in recorrelating some of the data reported 
for the three standard shapes”. In an attempt to improve this 
correlation, Krischer and Loos: © investigated the effect of 
fourteen shapes on the rate of heat and mass transfer. The 
results of their work will be discussed in more detail later; it 
should be noted at this time, however, that — because of the 
high value of the length-to-diameter ratio of the shapes they 
investigated — their characteristic length, L’, may be considered 
as useful for two-dimensional shapes only. 

In the present investigation three-dimensional shapes were 
studied with the following objectives: (1) To obtain convective 
heat and mass transfer data for stationary spheres and to com- 
pare the results with those in the literature; and (2) To deter- 
mine and compare convective heat and mass transfer data for 
other shapes in various orientations in a fluid stream in order 
to define a new characteristic dimension which will depend not 
only on the particle shape but on its orientation as well. 


EXPERIMENTAL 


The rates of heat and mass transfer during the evaporation 
of water from different stationary particles were studied in a 
| 1/2-in. diameter pyrex glass column in which turbulent air 
was passed under controlled conditions of temperature, humi- 
dity and velocity. 


Particle Shapes 


After much investigation, Celite was selected as the material 
from which the particles were made, because of its strength, 
inertness, and workability. The dimensions of the particles 
remained completely unaffected by moisture or air temper- 
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ature. Since Celite is composed of a microporous structure of 
minute silica particles — approximately four microns in size — 
held together by a bentonite clay binder, it loses most of its 
moisture at constant rate. An important added feature is that 
Celite changes color from orange to a pale yellow on drying, 
allowing a qualitativ e picture of the rate of mass transfer around 
a particular shape to be obtained. 

The various shapes investigated in this study, their dimen- 
sions and orientation in the flow stream, are given in Table 1. 


Apparatus 

Figure 1 shows a schematic diagram of the apparatus. It 
consisted essentially of two 1}4-in. i.d., 12-in. long Pyrex glass 
columns in parallel, through which hot air could be passed, in- 
stalled in a large electrically -heated, constant- temperature oil 
bath equipped with heaters, a cooling coil and a stirrer. The 
hot air from the }-in. supply line entered the bottom of one of 
the columns and passed in succession through a I-in. diameter 
orifice and a 14-mesh wire screen, 24 inches apart. A flat 


TABLE 


ALUMINUM 
=— ROO 


ene THERMOMETER 
~ COOLING 
COIL 
| \ HEATING 
COIL 


PYREX GLASS 
COLUMN 


——-OIL BATH 


6000 


WATT 
HEATER 


ROTAMETER 
THERMOMETER 
PRESSURE 

GAU 


GE \ 
©) 





PRESSURE - REDUCING VALVE 


Figure 1—Schematic diagram of apparatus for stationary 
shapes. 
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Panes. E SHAPES AND ) DIMENSIONS © 


Shape and Orientation 


Particle Designation and Dimensions 


Results 


Si: 
S2: 
S3: 


Cylinders 
(Disk) 


Cylinders 


Prism 


Hemispheres 


Figure 4 


0.380 in.; L 
0.197 in.; 
0.378 in.; 


Figure 5(A) 


0.264 in.; 
0.266 in.; 
0.264 in.; 
0.481 in.; 


Figure 5(B) 
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0.410 in. 


0.379 in.; 
0.150 in. 


Figure 6(A) 


0.360 in. 
0.252 in. 


Figure 6(A) 


0.282 in. 
0.244 in, 


Figure 6(B) 


0.284 in. Figure 6(B) 


0.247 in. 
0.251 in. 
0.222 in. 
0.247 in. 


Figure 7 
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velocity profile was thus obtained two inches from the top of 
the column, where each of the particles was suspended. 


A particle under test was suspended in place at the end of a 
3-in. long hypodermic needle held rigidly at the end of a small 
pow rerful permanent magnet embedded into the end of a long 
aluminum rod. In this way the particle remained free from vibra- 
tions in the air stream, and yet could be removed at will for 
weighing. 

The air was supplied by a large compressor, and passed in 
succession through a pressure-reducing valve, a calibrated rota- 
meter, a 6000-watt electrical heater and finally into a large heated 
surge tank from which it was led to the Pyrex columns through 
the insulated }-in. supply line. Its temperature was maintained 
at a preset value by an indicating temperature controller which 
regulated the electrical input to the heaters. In any one experi- 
mental run, the temperature of the air passing through the Pyrex 
column and the temperature of the oil bath were held at the same 
value to eliminate the occurrence of natural convection currents 
which would have been caused by cooler column walls. 


Procedure 


After completely steady conditions of air velocity and tem- 
perature had been attained in the apparatus, a Celite particle of 
known shape, weight and moisture content was suspended in a 
definite orientation in the glass column. As drying proceeded, 
the particle was weighed every few minutes by means of a deli- 
cate quartz-spiral balance, the extension of which was measured 
through a telescope. Each run was continued until the orange 
(wet) Celite turned a pale yellow color, indicating that the 
drying had passed from the constant-rate into the falling-rate 
period. The wet- and dry-bulb temperatures of the air in the 
column were accurately measured before and after each run so 
that the wet-bulb temperature of the particle as well as the air 
humidity could be recorded. 


Turbulence Intensity of Air Stream 


The velocity and turbulence intensity profiles for air passing 
through the empty column were obtained by means of a Model 
HWB2 Flow Corporation hot wire anemometer equipped with 
a HWP-B tungsten wire (0.00035-in. diameter) probe. The 
results obtained at two air velocities are shown in Table 2. The 
determinations were made at a cross-section 2}-in. below the 
top of the 12-in. long column and }-in. below the position where 
a particle was normally suspended in the column. The flatness 
of the velocity profile should be noted. 


TABLE 2 
AiR VELOCITY AND TURBULENCE INTENSITY PROFILES 


Fraction Air Intensity of Air Intensity of 
of column | velocity | turbulence | velocity | turbulence 
diameter ft./sec. Ye ft./sec. % 





1 10.6 10.3 
0.50 23.1 a 10.8 10.0 
0.75 22.4 _ 11.6 ~ 


0. 25 22.1 o. 
9 








Treatment of Experimental Data 


The results from each experimental run were processed so 
that they could be used in the form of a modified Chilton- 
Colburn equation “?); 
kcPiu 

Gm 


2p 


jo = (Sc)? = a(Re’”)-* 


Where a and b are constants, and Re’ is the particle Reynolds 
Number, based on its characteristic dimension. 

The reason for this choice of correlation was twofold: (a) 
it may be seen by rearranging the variables, that a plot of 
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the modified Nusselt Number for mass transfer, Nu’, as a func- 
tion of Re’’ is in reality a plot of (j ‘») (Re’’) (S89 against 
Re’’, so that it is more direct to plot j’s against Re’’ immed- 
iately and (b) the characteristic dimension, or diameter d, 
appears on both sides of the Nusselt equation while it appears 
only in the Reynolds Number in the Chilton-Colburn expres- 
sion. This makes the data for the latter equation easier to re- 
arrange and plot, especially in view of the fact that the object 
of this investigation was to find the value of the characteristic 
dimension. 

In determining Re’’, values. of viscosity and density were 
used at the average temperature of the gas film around the 
particle. Average values for the air velocity across the column 
were also used. For j’p, the mass transfer coefficient kg was 
calculated from the loss in water per unit time, the mass of 
the particle as well as the driving force involved between the 
vapor pressure of the water at the wet-bulb temperature at the 
surface of the particle and the partial pressure of water in the 
air stream. The Schmidt Number, Sc, was verified to be con- 
stant over the range of temperatures used; the diffusivity data 
for water evaporating in air being obtained from the Inter- 
national Critical Tables “®). 


Experimental Results 


One hundred and sixteen runs were carried out on the 20 
shapes described in Table 1 over the following range of operat- 
ing conditions: air velocities, from 4 to 25 ft./sec.; air temper- 
atures, from 70° to 250°F.; air humidities, from 0.003 to 0.010 
Ib. water/Ib. dry air. 


Figure 2 shows the qualitative observations obtained from 
changes in color of the Celite during the drying process of the 
various shapes. On this diagram the distance between the shape 
outline and the corresponding envelope at a particular point 
gives a measure of the relative rate of heat and mass transfer 
at that point. Shaded parts show areas of minimum transfer. 


Figure 3 shows the calculated results plotted (with the excep- 
tion of those for cylinders with their lengths parallel to the air 
stream) in the form of Equation (1) — with a possible error of 
+ 15% — using-a new particle characteristic dimension, L’’, 
instead of the usual diameter, d, in the Reynolds Number. This 
new characteristic dimension, L’’, for a stationary particle of 
any shape suspended in a moving fluid stream in a definite orien- 
tation is defined as the total surface area of the particle divided 
by the perimeter of the maximum projected area perpendicular 
to flow. As an example, for a cylinder of length, L, and dia- 
meter, d, suspended with its length perpendicular to flow: 


Lr surface area 
max. perimeter 


wdadL + 2md?/4 
2(L + d) 


When the cylinder has a high aspect ratio, i.e., when L is much 
greater than d: 
ad 


L” = — = L’ 
2 


where L’ is the characteristic “flow length’ as defined by 
Krischer and Loos. 


Discussion of Results 
Spheres 
Figure 4 shows the results for four spheres plotted snd cor- 


related by the following generalized equation with an average 
deviation of + 10% (for a sphere, L’’ = d and Re’’ = Re): 


= 0.692 (Re’’)-?.488 


Since no special trend in the scatter of points can be noted 
for each of the comparatively large spheres used, the wall effect 
on the rate of heat and mass transfer was considered to be 
negligible or constant — within the experimental deviation. 
Because the size range for all other shapes was similar to that of 
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Figure 2 — Qualitative transfer rates around the shapes 
used. 


-0.486 


a 0.692 (RE") 


600 1000 2000 4000 000 
Re" 


Figure 3—Results for shapes used. 


400 600 1000 2000 4000 10,000 20,000 40,000 60,000 
RE 


Figure 4—Data for spheres. 


spheres, the wall effect was neglected throughout the 
investigation. 

Figure 4 also shows the data of Powell, Williams“, Maisel 
and Sherwood“, Linton and Sherwood“® and Krischer and 
Loos®: , It may be seen that the proposed correlation differs 
slightly in slope from the experimental lines of the other invest- 
igators but falls between them. The scatter in the data from 
the other investigators may be partly due to the fact that in 
some cases the average velocity across the column was used — 
as in this investigation — while in other cases only the approach 
velocity to the particle was measured; it may also be partly 
due to the fact that each investigation was carried out under 
different conditions of fluid turbulence. Only recently have 
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attempts been made to correlate the turbulence parameters with 
heat and mass transfer rates “5: 2 2, 22), 


Qualitatively it was observed that the maximum rate of 


heat and mass transfer occurred over the front half of the sphere 
with a second maximum occurring at the rear — around the 
point of suspension. A wet ring around the sphere at from 
100 to 110° from the stagnation point indicated the position of 
minimum transfer. These observations agree with those of 
Cary) as well as those of Garner) who noted that for a solid 
sphere at Reynolds Numbers greater than 500 the minimum 
mass transfer rate occurred at approximately 105° from the 
stagnation point. 


Cylinders 


(a) Length perpendicular to flow. Figure 5 (A) shows 
that not only have cylinders of different length-to-diameter ratios 
been correlated with spheres, but a disk (C3) also appears to 
follow the same behavior. 

The heat and mass transfer flux around a cylinder in this 
orientation, as shown in Figure 2, agreed well with the observ- 
ations of Robinson and Han®® and Winding and Cheney °® who 
noted a maximum Nusselt Number at the front and rear of the 
cylinder with a minimum at approximately 90° from the stag- 
nation point. The rate of heat and mass transfer from the two 
flat sides of the cylinder was intermediate between the minimum 
and maximum values. 


(b) Length parallel to flow. Figure 5(B) shows that the 
correlation line does not follow the line for spheres completely 
but breaks away at Re’’ = 1000, presumably because at higher 
values the boundary layer becomes turbulent. This explanation 
is based on a possible analogy with the plot given in the Chemical 
Engineers’ Handbook” for the rate of heat and mass transfer 
from plane and cylindrical surfaces located parallel to the 
direction of flow. This plot compared the data of Pasquill @®, 
Wade) and Powell with the theoretical equations of 
Pohlhausen™ and Latzko®, and clearly establishes the regions 
where the laminar and turbulent transfer rates occur. In this 
investigation, early transition to a turbulent boundary layer 
may well be caused by the sharp edge around the upstream 
diameter. 

It was observed that for a cylinder in this orientation the 
position of the minimum transfer point, which occurred on a 
circular band around the rear portion (see Figure 2), did vary 
with air velocity. However, no definite conclusions could be 
drawn regarding its changes in location with increasing Reynolds 
Numbers. An increase in the relative heat and mass transfer 
rate around the edges is probably due to a thinning of the 
boundary layer in these regions. 


Prism 


Only six experimental points were obtained for a flat prism 
suspended with its faces parallel to the air stream and the results 
shown in Figure 6(A) are inconclusive. It cannot be ascertained 
whether the single point below the sphere correlation is due to 
an experimental error while the first five points indicate a 
tendency for a turbulent boundary layer to be formed, or 
whether all the points just show an experimental scatter of + 16% 
about the sphere line. The latter case is also possible since it 
was extremely difficult to suspend a prism with its faces exactly 
parallel to the air stream. Qualitatively the maximum transfer 
rates occurred at the front with the minimum at the rear, as 
shown in Figure 2. Again high transfer around the edges was 
noted. 


Cubes 


(a) Face parallel to air stream. |n this case Figure 6(.4) 
shows that all results follow the correlation for spheres quite 
faithfully in the range of Reynolds Numbers and for the two 
cube sizes investigated. 
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(b) Corner facing the air stream. The results in this case 
follow the line for spheres but are somewhat lower, as shown 
in Figure 6(B). The reason could well reside in the suspension 
method, as any departure from the exact position would cause a 
decrease in the actual value of the perimeter, thus making the 
value of L’’ larger than that calculated. 


(c) Edge facing the air stream. Figure 6(B) shows that 
the results are entirely similar to the previous ones, quite possibly 
for the same reason. 


vr 


For the cube in case (a) : L’’ = 1. 

For the cube in case (b) : L’’ = 1. 

For the cube in case (c) : L’’ = 1. 
where L is the length of one side of the cube. 
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The transfer fluxes for the cubes in the three orientations 
were observed to be approximately as shown in Figure 2. For 
all three cubes the rate of heat and mass transfer around the 
edges appeared to be very high relative to the other areas. The 
minimum point of transfer is noted by the shaded area. 


Hemisphere 

Figure 7 clearly establishes the validity of using L’’ for 
hemispheres (H1 and H2), or part of a hemisphere (H3). 
However, for hemisphere Hla, in which the flat rear section 
was covered with an impervious layer of Pliabond cement, 
the data fall about 40-50% higher than the sphere line. This 
behavior was to be expected in view of the fact that previous 
published results for a sphere have shown that the ratio of the 
mass transfer before the point of separation to that after the 
point of separation is of the order of 75 to 25, depending on 
the magnitude of the Reynolds Number®). From Figure 7 
it may be shown that the ratio of mass transfer on the curved 
surface of the hemisphere (before the point of separation) to 
that on the flat surface (after the point of separation) is 96.5 
to 3.5 at Re’’ = 1000 and 87.1 to 12.9 at Re’’ = 3000. These 
ratios are of course independent of the characteristic dimension 
used for the particle diameter. 


Correlation of Data in Literature Using L” 


Further confirmation for the validity of Equation (3) was 
sought from the published literature. A search revealed few 
investigations where more than one shape had been used. The 
six major ones which were thus found were correlated and 
replotted in the form of j’p — Re’’ graphs. These graphs are 


shown in Figures 8 to 11 inclusive. 


Data of Maisel and Sherwood(!*) 


Maisel and Sherwood reported their results for spheres, 
cylinders (axis perpendicular to air flow) and disks (diameter 
perpendicular to air flow) in the form of three separate j’» — Re 
graphs. Since they were mainly interested in obtaining accurate 
evaporation data for various volatile liquids, no attempt was 
made to present a unified correlation. Figure 8(A) shows their 
data correlated with that of a sphere using L’’ as the character- 
istic dimension. It may be seen that the agreement is excellent 
and that a single correlation could replace their three plots 
quite effectively. 


Data of Linton and Sherwood(!) 


Linton and Sherwood also reported their results for spheres, 
cylinders and flat plates in the form of separate j’p— Re graphs. 
Figure 8(B) shows their results for the benzoic acid sphere and 
cylinder and the cinnamic acid flat plates correlated by the use 
of L’’. The change in slope for the flat plates, indicating the 
formation of a turbulent boundary layer, should be noted. 


Data of Norris and Spofford(32) 


These authors presented in 1942 a generalized correlation 
for convective heat transfer from pin fins and strip fips of various 
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designs, and compared their results with the data of Fage and 
Falkner for a single plane, the data of Colburn“ for a 
single cylinder with its length perpendicular to the air flow, 
and the data of Grimison®® for a bank of cylinders. They 
presented their correlation and comparison data in the form of a 
ju— Re, graph and obtained the equation: 


ja = (Re, f°... 555s 


where jq—heat transfer factor; 

Re,—Reynolds Number based on the perimeter, P, of the 
shapes. 

It may be shown that for all the shapes used and correlated, 
the authors were using P = 2L’’ (approximately) as their 
shape parameter. Figure 9 shows their graph replotted using 
L’’ in the Reynolds Number. The agreement in correlation 
and comparison is good. 


Data of Hilpert(?®) 

Figure 10 shows Hilpert’s data for a long cylinder and long 
tubes of various profiles, replotted in terms of L’’ and compared 
with Colburn’s equation for a flat plate. Again the correlation 
appears to be good. This graph indicates the variations that 
may be expected when dealing with bodies of different shape. 


Data of Powell(®) 


Figure » shows Powell’s data replotted and correlated by 
means of L’’, assuming an air temperature of 20°C. The heat 
and mass transfer coefficients for the circular disks and rec- 
tangular plates facing the air stream (lines E-E and G-G) were 
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obtained with water evaporating from the front surface only. 
However, in calculating L’’ for these shapes, the total surface 
area of both sides of the thin shapes was used rather than just 
the “wetted” area. The resulting correlation shows that the 
type of flow around a particle is dependent only on the particle 
dimensions (as it should be) and is independent of the area 
active for heat and mass transfer. When the wetted areas 
faced downstream, the heat and mass transfer rates were from 
10-40% lower. 


Data of Krischer and Loos(®,1°) 


In 1958 Krischer and Loos obtained a definition for the 
characteristic length, L’, of an irregular particle by studying 
the rate of heat and mass transfer from 14 shapes. All these 
shapes were long and thin, 1.e., they had a high aspect ratio. 
On the basis of this rather restricted particle geometry, they 
defined L’ as the length of one side of a rectangle whose surface 
area is equal to one half the surface area of the body and whose 
other side has a length equal to the length of the shape perpendic- 
ular to the direction of flow. For a shape whose surfaces were 
only partially involved in the transfer processes, the average L' 
was taken as half the length of the development of the exchange 
area. 


It should be noted that for all the shapes used by Krischer and 
Loos except spheres, L’’ is approximately equal to L’. However, 
L"’ for a particular shape appears to be less arbitrary, easier 
to define and more practical to use. More important still, 
L’ does not apply to a particle shape whose dimensions along 
all three axes are roughly similar. ‘Thus for a sphere L’ = d/2, 
while L’’ = d. 

Krischer and Loos’ following general conclusions are worthy 
of note: (a) L’ can be used as the characteristic length for any 
given shape (with a high aspect ratio) with sufficient exactness 
to correlate the heat and mass transfer coefficients with the 
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Reynolds Number (+ 15 percent) provided that no recessed 
surfaces are present. (b) For bodies with recessed surfaces, 
low heat and mass transfer results are obtained at low Reynolds 
Numbers (<2500) but approach the values for standard shapes 
(spheres, cylinders and flat plates) at high Reynolds Numbers 
(>2500). The low values are due to dead areas of fluid with 
little turbulence in the cavities. (c) For a body whose true 
surface and flow length L’ are difficult to define, the “sheath 
area” and its flow length are better relative values to use than 
the surface area and flow length of a sphere or a cylinder of the 
same volume. (d) For all bodies with recessed surfaces, the 
sheath area appears to be the most suitable value for small 
Reynolds Numbers whereas for large Reynolds Numbers, the 
true surface area appears to be the better relative value for the 
determination of heat and mass transfer. 


Fundamental Significance of L” 


The present stage of development of the Boundary Layer 
Theory does not allow the derivation of a mathematical 
expression to show why L’’ should account for the three- 
dimensional effects of particle shape and orientation. It does, 
however, enable a qualitative understanding of its meaning to 
be obtained. 

From the general properties of a laminar boundary layer, it is 
known that the /ocal heat transfer coefficient h in the equation: 


dq = hdA (T, — T;) 


is inversely proportional to the thickness, 6;, of the thermal 
boundary layer, which is itself related to the thickness, 6,, of the 
velocity boundary layer and is nearly equal to it for a Prandtl 
Number of one. In addition, it is known that the Nusselt 
Number is proportional to (Re)!? and hence h (or kg) is in- 
versely proportional to the square root of the characteristic 
dimension used to define the Reynolds Number. However, to 
integrate Equation (5) for any body up to the point of separation 
(assuming constant temperatures) : 


the manner in which 6; varies as integration along the surface 
proceeds must also be known. This can be calculated for flat 
plates, regular and elliptical cylinders, spheres and a few other 
shapes 37), but the calculations become hopelessly complex for 
the shapes used in this investigation due to such factors as the 
possible separation and later re-attachment of the boundary 
layer at sharp forward edges. 

Since for a sphere 5; is anywhere proportional to d'? for 
a laminar boundary layer up to the point of separation (and 
proportional to d°-?° for a turbulent boundary layer) it is natural 
to take the diameter d of the sphere as its basic characteristic 
dimension for use in heat (and mass) transfer correlations. 
Once these correlations have been established for spherical 
particles (as was done in this investigation) it is logical to apply 
them in their original form to other “shapes by using the proper 
characteristic dimension. The problem thus resolves itself 
in finding, for any shape in any orientation, a characteristic 
dimension which — in an overall or integrated way (as distinct 
from local) — will account for the average thickness of the 
boundary layer. 

it can be verified — from the limited evidence available — 
that L’’ as defined appears to fulfill this requirement from both 
perimeter and surface area considerations. Thus, for shapes 
with a constant maximum perimeter, L’’ for any shape will 
increase as the surface area is increased, from its definition. 
But as the latter is increased, the boundary layer will thicken’, 
Figure 7.7) to cause an overall decrease in the integrated value 
of h (or kg). This is in agreement with the theory just presented 
that /) is inversely proportional to the square root of a character- 
istic dimension. From Equation (3) — obtained experimentally 
—it is evident that: 
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kg (overall) « 1/(L’’)-48* 


For shapes with a constant surface area, L’’ for any shape 
will increase as the perimeter is decreased, from its definition. 
However, when the latter is decreased, the boundary layer also 
thickens(*”, Figure 12.11) and an overall decrease in the inte- 
grated value of h (or kg) will result. Again agreement of 
experimental values with theory is obtained when L’’ is used 
as the characteristic dimension. 


CONCLUSIONS 


A new dimension, L’’, characteristic for single, stationary 
particles in a fluid stream, was introduced to account for both 
particle shape and orientation in the forced convection heat 
and mass transfer equations. This dimension is defined as the 
total surface area of the particle divided by the perimeter of the 
maximum projected area perpendicular to flow. 


2. The heat and mass transfer data for 20 shapes in different 
orientations were correlated with the aid of L’’ by the equations: 


J'p = 0.692 (Re’’)-°-486 for 500 < Re” < 5000 


with a deviation of + 15%. This equation can be readily 
expressed in terms of the Nusselt Number for mass transfer, 
Nu’, since Nu’ = j’p (Re’’) (Sc)"%, yielding the form: 
Nu’ = 0.692 (Re’’)%-5!4 (Sc)!/3, Zz can similarly be expressed 
in terms of the Nusselt Number, / 


3. The results for spheres were saad by this equation 
with a deviation of + 10%, and agreed quite w ell with published 
data. 


4. All the heat and mass transfer data in this investigation 
were obtained at an intensity of turbulence between 9 and 10% 


5. Qualitative pictures for the relative rates of heat and 
mass transfer around the various shapes used were presented 
(Figure 2). 

6. L’’ appears to be a more general shape parameter than 
Krischer and Loos’ L’. It was found to correlate the various 
data for shapes available in the literature — including those of 
Krischer and Loos. 


7. The general conclusions of Krischer and Loos concerning 
shapes with recessed surfaces or surfaces which are difficult 
to define probably apply when using L’’ instead of L’. but this 
assumption still requires experimental confirmation. 


8. A qualitative explanation for the meaning of L’’ was 
presented from boundary layer theory considerations. These 
considerations show that the use of L’’ as a characteristic 
dimension does account for the thickness of the boundary 
layer around a specific particle. 
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Nomenclature 


Roman Signals 


a and b = constants in Equation (1). 

/ surface area of particle, ft.? 
surface area of a sphere with same volume as particle, 
ft.? 
specific heat of air at constant pressure, B.t.u./(Ib.) 
(°F.). 
diameter, ft.; (d,, equivalent diameter). 
diameter of flat (back) part of hemisphere, ft. 
diffusion coefficient for water evaporating in air, 
ft.2/hr. 
mass velocity of flowing air, Ib./(ft.*)(hr.). 
molal mass velocity of flowing air, lb. moles/(ft.?)(hr.). 
heat transfer coefficient, B.t.u./(hr.)(ft.2)(°F.). 
average thermal conductivity of air film, B.t.u./(hr.) 
(ft.)(°F.). 
mass transfer coefficient, Ib. moles/(hr.)(ft.2)( Ap). 
constant in equation (6). 
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Ts 
r, 


length, ft. 

flow length of Krischer, ft. 

characteristic length of Krischer and Loos, ft. 

new characteristic dimension, ft. 

perimeter of a shape, ft. 

log mean pressure difference of air across the boundary 
layer, atmosphere. 

heat transferred, B.t.u./hr. 

average humid heat of the gas film, (B.t.u.)/(Ib. dry 
air)(°F.). 

arithmetic mean air film temperature, °F. 

dry bulb air temperature, °F. 

absolute temperature of bulk air flow, °R. 

absolute temperature of particle surface, °R. 


annua 


tnu ol 


Greek Symbols 


thickness of velocity boundary layer, ft. 
thickness of thermal boundary layer, ft. 
average viscosity of air film, lb. mass/(ft.)(hr.). 
average density of air film, lb./ft.* 

Leva’s shape factor, dimensionless. 

Carman’s shape factor, dimensionless. 

Wadell's sphericity = A,/A, dimensionless. 


Dimensionless Groups 


(2 


~ 


(3) 
(4) 


(5) 
(6) 


(7) 


(8 


~ 


(9 


ww 


(10) 
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Ip = Mass transfer factor = (kgPzy/G,,)(Sc)?!* 
ja = Heat transfer factor = (h/c,G)(Pr)?!* 
Nu = Nusselt Number = hd/k, (for spheres) or hL’’/k, 
Nu’ = Nusselt Number for mass transfer = 
kgML” Pry /D,py = Nu (Ry/sD,p,) 
Pr = Prandtl Number for heat transfer = c,uy;/ky 
Re = Reynolds Number for spheres = dG/py 
Re’ = Reynolds Number = L"’G/p, 
Re, = Reynolds Number = PG/wy; (equation 4) 
Sc = Schmidt Number for mass transfer = y,/D,p; 
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The Effect of Temperature on the 
Fluidization of Silver Powders’ 


E. ECHIGOYA? and G. L. OSBERG* 


Fluidization tests on pure silver and silver catalyst 
powder show that agglomeration occurs when the bed 
is heated, which finally results in bed collapse at a 
comparatively low temperature. If the bed temper- 
ature is held constant in the agglomerating temper- 
ature range, an induction period is observed before 
fluidization ceases. A study of the induction period 
as a function of temperature (10-60°C.) leads to an 
activation energy of 18 keal. for the overall agglo- 
meration process for pure silver. Agglomeration is 
believed to be due to sintering. 


NTEREsT in these laboratories in fluidized silver powders 

developed from exploratory design studies on catalysts 
prepared from Ca Ag alloys “) which have been proposed 
for the oxidation of ethylene to ethylene oxide. Small 
scale fluidizing tests with “these catalysts showed that the 
formation of agglomerates, when the ‘bed temperature was 
raised, prevented stable fluidization in the operating 
temperature range. When fluidization ceased, catalyst 
tempeartures were uncontrollable because of poor heat 
transfer. Since little information has appeared on the 
agglomeration of fluidized metal powders, it is perhaps 
of interest to summarize our observations on fluidized 
silver powder systems. 


Experimental Method 


The fluidizing tests were carried out in a 2.5 cm. LD. 
or 5 cm. L.D., electrically heated, pyrex tube with a 
conical (23° included angie) bottom through which the 
fluidizing gas was fed. The inlet gas passed through a 
drying tube and rotameter, and the outlet gas through a 
cyclone. A manometer was connected to the bottom of 
the cone, and to the top of the pyrex tube, so that 
pressure drop measurement could be taken across the bed. 
Bed temperatures were measured with a thermocouple 
dipping into the bed. 

Fluidizing difficulties due to agglomerate formation 
were directly visible through the pyrex tube, but were 
also readily deduced from pressure drop readings. With 
normal fluidization, pressure drop readings fluctuated 
slightly, but with the mean value equal to the weight 
of the bed. As agglomeration occurred, parts of the bed 
were immobilized; channels developed through the bed, 
causing a sudden decrease in pressure drop. Under some 
conditions these large erratic fluctuations were repeated 
several times before the bed was completely immobilized. 
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2N.R.C. Postdoctorate Fellow. Present address: Tokyo Institute of Tech- 
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Powder samples used in these experiments were silver 
powder and silver catalyst. The silver powder was of 
high purity containing .0005% Pb, .005% SiO, .0004% Fe, 
.0002% Al. The silver catalyst was prepared from 8.5% 
Ca, 91.5% Ag alloy powder as described by McKim and 
Cambron “), The calcium was removed by leaching with 
20% acetic acid after pre-treatment for one hour with 
water in an autoclave at 225°C. Following leaching, wash- 
ing and drying, a catalytically-active silver catalyst was 
obtained having a residual calcium content of about 0.15% 
and a specific surface area of 0.5 sq. m./gm. 


Observations 


At room temperature silver catalyst powders fluidized 
satisfactorily in air, but when heated ceased to fluidize at 
a definite temperature. In Figure 1, the arrest temperatures 
are given for three mesh sizes (Tyler) of silver catalyst. 
It may be noted that the finer powders show substantially 
lower arrest temperatures, and that air flow rates within 
the range shown had only a minor effect. A separate 
series of tests with sample sizes varying from 10-25 ce. in 
the 2.5 cm. LD. tube, showed no significant effect of 
volume of sample on the arrest temperature observed. 

The agglomerates which formed at these low temper- 
atures were of relatively low strength. Cooling the bed to 
room temperature, and maintaining the gas flow was often 
sufficient to bring about fluidization again. If this cycle of 
heating and cooling was repeated several times, however, 
a lower temperature was observed after each cycle 
though the arrest temperature appeared to be tending 
toward a constant value. 


Induction Period 


When a fresh sample of powder was maintained at a 
constant temperature under fluidizing conditions diffi- 
culties due to agglomeration did not occur immediately. 
A time interval was observed which depended on the bed 
temperature. An induction period was therefore selected as 
the time interval between the beginning of fluidization, 
and the first major decrease in pressure drop which 
culminated in bed immobility. Induction periods, thus 
defined, were observed with pure silver, and with silver 
catalyst powders using air, hydrogen and argon as fluidiz- 
ing gases. When hydrogen was used there was no difficulty 
in observing the induction period because bed immobility 
occurred immediately after the first major decrease in 
pressure. When air or argon was used, major pressure 
fluctuations were erratic and occurred several times, with 
normal readings in between, before fluidization finally 
ceased. For this reason, the first major fluctuation in the 
series which culminated in immobilization of the bed was 
arbitrarily selected as marking the end of the induction 
period. 
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Figure 1 — Fluidizing arrest temperatures for different 
mesh ranges of silver catalyst in air in a 2.5 cm. I.D. 
fluidizer. 


The quality of fluidization also had an effect on the 
observed induction period. When slugging conditions 
were encountered, longer induction periods were observed. 
With small bed volumes, i.e. 5 to 12 cc., slugging was 
not encountered, and reproducible induction periods were 
obtained. 

Induction period data are given in Figure 2 for a 
temperature range of 10 to 60°C. These bed temperatures 
were obtained by immersing the glass fluidizer in a 
constant temperature water bath. Bed sample size (10 cc.) 
and gas flow rates were constant for a particular set of 
temperature data. Gas flow rates at N.T.P. to the 2.5 cm. 
1.D. fluidizer were 2.8 litre/min. for air and argon, and 
3.7 litre/min. for hydrogen. 

The Arrhenius plot satisfactorily correlated the induc- 
tion period as a function of reciprocal bed temperature 
so that the data may be represented by an empirical 
equation 

1 = A e@~AE/RT 

t 
where t¢ is the induction period, min., T is the bed tem- 
perature, °K., and A and E are experimental constants. 
Using the usual Arrhenius interpretation, AE is regarded 
as the activation energy for the overall agglomeration 


process. 
Two groups of approximately parallel lines are shown 


in Figure 2. Plots A, B, and C were obtained with pure 
silver as the bed material. These have approximately the 
same slope, but have different intercepts. Changes in par- 
ticle size, fluidizing gas and history of pure silver caused 
variations in the intercept, i.e. constant A. Changes in 
quality of fluidization which were obtained by using a 
5 cm. I.D. glass unit instead of the 2.5 cm. I.D. unit 
reported here, also caused a change in the value of the 
intercept but not in the slope. The slope corresponds to 
an activation energy of about 18 kcal. 

The second group of lines, D. E., and F, were obtained 
with silver catalyst as the bed material. Lines D and E 
are almost parallel while F is somewhat steeper. All are 


Ad 
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Figure 2—Induction period ys reciprocal temperature. 


at a steeper slope than those for pure silver. The effect 
of mesh size, etc., is as above reflected in the intercept. 
The activation energy calculated from slope of these 
lines varies from 38 to 52 kcal. 


Discussion of Data 

It is evident from these limited data that the overall 
agglomeration process is a complex one. Fluidizing vari- 
ables such as gas velocity, particle size, column design, 
etc., as well as surface properties of the metal powders 
all hove an effect on the observed agglomeration rate, and 
no attempt is made to develop general correlations. 

The formation of agglomerates is believed to be due 
to low temperature sintering of particles when these 
particles come into contact following collision within the 
fluidized bed. Low temperature sintering of metal powders 
is well known, and has been studied by Tammann and 
Mansuri) in a similar way, by noting the temperature 
at which stirrer ceases to move when stirring a metal 
powder which is being heated. They report arrest tem- 
peratures for silver in the range of 136-145°C. Thus, the 
fluidized bed tests above may be regarded as a different 
technique, perhaps a more sensitive one, for establishing 
the onset of sintering. 

The mechanism of sintering of metal powder under 
low stress has received limited study. Kuczynski () has 
suggested that there is a close connection between sinter- 
ing and the process of self diffusion. He observed experi- 
mentally that the rate of interfacial growth between a 
copper sphere in contact with a copper sheet is deter- 
mined by the volume diffusion rate at high temperatures 
and predominantly by a surface diffusion mechanism at 
a low temperature. 

Following this suggestion, it is interesting to compare 
the observed activation energy of 18 kcal. for agglomer- 
ating pure silver powder with published values of diffusion 
processes in silver. Hoffman and Turnbull () report an 
activation energy of 20.2 kcal. for grain boundary data 
obtained on polycrystalline silver at 400°C. They give a 
value of 45.9 kcal. for the lattice diffusion rate in a single 
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crystal of silver. Nickerson and Parker (°) report an acti- 
vation energy of 10.3 kcal. for the surface diffusion of 
silver. [he activation energy obtained from the agglome- 
ration data is intermediate between the surface and lattice 
diffusion values, and rather close to the grain boundary 
diffusion value. The activation energy of 38-50 kcal. for 
silver catalyst would appear to be higher than for pure 


by the addition of small amounts of CaCO, “), and by 
modification of the basic cal: um silver alloy through the 
addition of hardening elements such as nickel. All these 
methods have been effective to a degree in overcoming 
agglomeration, but none, so far, have been completely 
satisfactory from both the fluidizing, and the catalytic 
point of view. 


silver because of the presence of the impurity CaO, which 
might be regarded as hindering the mobility of the References 
diffusing silver. The data presented here are, however, (1) McKim, F., Cambron, A., U.S. Pat. 2,562,857 (1951). Can. 


too limited to justify detailed speculations on the mechan- Pat. 475,366 (1951). 

ism of sintering. Tammann, G., Mansuri, O. A., Anorg. Allgem. Chemie 
126, 119 (1923). 

Kuczynski, G. C., Journal of Metals 1, 169 (1949). 

Hoffman, R. E., Turnbull, D., J. Appl. Phys. 22, 634 

(1951). 

Nickerson, R. A., Parker, E. R., Trans. Am. Soc. Metals 

42, 376 (1950). 


Minimizing Agglomerate Formation 

These data show that the effect of agglomerate forma- 
tion on fluidization can be minimized if coarse rather than 
fine silver particles are used, and if surface impurities such 
as CaO are present. Attempts have been made to render Sears, G. W., U.S. Pat. 2,615,899. 
the silver catalysts non-agglomerating by the addition of Osberg, G. L., U.S. Pat. 2,537,486. 
diluent powders such as graphite as proposed by Sears (°), -- * 


Edmonton, Alta. 


Editor, C.J.Ch.E.: 


We wish to point out the following errors which 
appeared in the February 1959 issue of the C.J.Ch.E. We 
thank B. C. Sakiadis of E.I. Du Pont de Nemours and 


Company for bringing them to our attention. 


1. In the article “Horizontal Pipeline Flow of Mix- 
tures of Oil and Water” by T. W. F. Russell, G. W. 
Hodgson and G. W. Govier, the term 8ay in the defini- 
tion of M, in the Nomenclature should read 8ay3. 

2. In the article “The Effect of the Less Viscous 
Liquid in the Laminar Flow of Two Immiscible Liquids” 
by T. W. F. Russell and M. E. Charles the term 8ay in 
the definition of M, in Equation (7) should read 8ay3 
and the terms (3u, + 6u, + Ilu,ug) a*y® in the equa- 
tion between Equations (9) and (10) should read 4(3u,? 
+ 6ui, — 11uaupg)a*y?. 


G. W. Hodgson and M. E. Charles, 
Research Council of Alberta. 
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A Method for the Determination of 
Tie-Lines in Ternary Liquid Systems: 


E. L. HERIC? and RONALD M. RUTLEDGE: 


A method is described for the determination of 
tie-lines in ternary liquid systems which may be used 
when a satisfactory method of chemical analysis is 
unavailable. The technique is graphical and is based 
upon the lever rule principle. The method is illus- 
trated with the system acetic acid-furfural-water at 
25 + 0.05°. Results indicate the method is accurate 
and capable of good precision. 


r a study of the system acetic acid-furfural-water in a manner 
rather more refined than those used previously “:*), a suitable 
method was needed for the determination of the tie-lines. Usual- 
ly the presence of an acid as one of the components suggests 
titration of that component with a base as the simplest method. 
In the present instance, however, it was felt that this method 
might not be feasible because of the instability of furfural ®*). 
It was thought that the acidic nature of some of the decompo- 
sition products might introduce a systematic error if the method 
of tie-line determination were based upon acid titration. 

There are available several methods of tie-line determination 
based upon the lever rule of Miller and McPherson“. The 
methods of Evans‘? and Othmer and Tobias‘ both require a 
knowledge of the total weights of a synthetic sample and each 
of the conjugate phases formed from the sample. Another tech- 
nique described by Bancroft and Hubard‘* involves working 
with only a representative fraction of each of the conjugate 
phases. Johnson and Barry“), using a mechanical device, mod- 
ified the latter method by eliminating some of its undesirable 
time-consuming trial and error features. While the approach of 
Johnson and Barry is convenient, it was hoped that our data 
would be of sufficient precision to warrant the use of a larger 
graphical scale for evaluation, rendering the size of their mechan- 
ical device impractically large. The method described in this 

aper is another approach to modification of the method of 
Bancroft and Hubard, again removing its trial and error features 
but requiring no mechanical device other than a straight-edge. 


Chemicals and Equipment 


Furfural, commercial grade, was puritied by distillation under 
vacuum, and only the middle fractions were used. For freshly 
distilled furfural, mp at 25 + 0.05° was 1.5235. The accepted 
value is 1.523459. To minimize decomposition of the furfural 
during storage, it was kept in closed containers at the sublimation 
temperature of carbon dioxide. When used, each sample of 
furfural was colorless, although a slight coloration was frequent- 


1Manuscript received August 19, 1959; Accepted January 15, 1960. 
2Assistant Professor, Department of Chemistry, University of Georgia, 
Athens, Georgia. 


3Present address: Department of Chemistry, University of Tennessee, 
Knoxville, Tennessee. 
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ly evident by the time the sample was discarded. Acetic acid 
was of reagent grade quality, minimum assay 99.7% by weight. 
It was used without further purification. Water was from a 
laboratory still. The thermostated samples were maintained at a 


temperature of 25 + 0.05°. 


Determination of the Binodal Curve 

The binodal curve was determined by adding acetic acid to 
synthetic samples in the two-phase region until homogeneity 
was established. From the weights of the three components 
present, points on the curve were established. Because of the 
instability of furfural only one or, at the most, two points on 
the curve were established with a given sample of furfural. 
While duplicate determinations were not made for three-com- 
ponent samples, an indication of the precision of this work is 
given by the results obtained in repeated determinations of the 
solubilities of furfural in water and water in furfural. In both 
the former, where three determinations were made, and the 
latter, where four determinations were made, the average devia- 
tions were less than 0.1% by weight. The solubility curve of 
the system is given in Figure 1. 

In passing, a remark concerning the effect of furfural decom- 
position on solubility is pertinent here. It has been found that 
the decomposition of furfural has little effect on solubility in 
some instances“, However, when compared with the present 
work, data obtained by one of the authors previously”) on the 
solubility of furfural in water and vice versa (with less carefully 
purified furfural than used here) shows sufficient disagreement 
with the present work to throw doubt upon the generality of 
such an assumption. 


Determination of the Tie-Lines 


Tie-line determination was accomplished with the aid of 
Figures 2 and 3. The former is a reconstruction of Figure 1 in 
terms of a right triangle*, and is necessary because the scale of 
the usual equilateral triangular graph paper is too small for 
precise work. By adopting a right triangle it is possible to use 
an extended form of rectangular coordinate paper, giving a 
higher degree of sensitivity to the method. (Figure 2 has been 
drawn out of scale in order to show more clearly the essential 
features of the method.) 

Points K and L in Figure 2 represent two binary mixtures 
of respectively acetic acid and furfural, and acetic acid and water. 
From K a number of straight lines were drawn through the 
binodal curve, such as KN,Z1, KN2Z. and KN3Z;3. The loca- 
tions of the points such as Z,, Z2, Z; are arbitrary, except that 


*In practice the construction of Figure 1 is not necessary for the applica- 
tion of the method to be described. It is used here only as a means of 
presenting the data of the system in the more usually encountered form. 
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in the present system, for example, they were to cover the 
range i solutions found on the water-rich portion of the binodal 
curve. Ina like manner straight lines were drawn through the 
binodai curve from point L, such as YM ,L, ¥2M2L and Y3M3L. 
Again che locations of the points such as Y;, ¥2, Ys are arbitrary, 
except that in the present system they were to cover the range 
of solutions found on the furfural-rich portion of the binodal 
curve. Because of the arbitrary way in which the Y and Z 

ints have been selected, it is apparent that like subscripts on a 
pair of points Y and Z do not infer that these points bear any 
relationship to one another. At each of the points N and Z on 
a given line KNZ the weight percentages of water, Wy and Wz, 
were read from Figure 2 and recorded. It was thus possible to 
construct the curve in Figure 3 of Wy/Wz, versus Wz. Similarly 
at each of the points Y and M on a given line YML the weight 
percentages of furfural, Fy and Fy, were read from Figure 2 
and recorded. The curve in Figure 3 of Fy/Fy versus Fy was 
then constructed. By combining the information in Figure 3 
with the experimental procedure described next, it was possible 
to determine the tie-lines of the system. 

A series of synthetic samples were prepared. Each sample 
separated into two conjugate phases, Y and Z, whose composi- 
tions were to be determined. A fraction of each of the two 
phases in each sample was removed and titrated with a titrant 
solution. The composition of the titrant solution corresponded 
to either K or L in Figure 2. Solution K, containing 20.0 weight 
% acetic acid and 80.0 weight % furfural, was used to titrate 
those fractions on the water rich side of the binodal curve, the Z 
phases. Solution L, containing 20.0 weight % acetic acid and 
80.0 weight % water, was used to titrate those fractions on the 
furfural rich side of the binodal curve, the Y phases. In every 
case the fractions were titrated to clearness following the initial ACETIC 
formation of a cloudy suspension. At the “end points’’ of the 
titrations the compositions of the titrated fractions were points 
M and N on the opposite side of the binodal curve from the 
fraction being titrated. M and N lie on the straight lines joining 
L and ¥ in the one instance, and K and Z in the other. 

By the lever rule, if to a weight gz of composition Z there 
is added a weight gx of composition K in order to achieve a 
composition of N, then 


Figure 1—The system acetic acid-furfural-water at 25°. 


NZ 


Figure 2—An illustration of the method used in construct- 


Moreover it follows from Equation (1) that : : da 
‘ q (1) ing the working plot for tie-line determination. 


» _&R 
&x + gz ok 


In addition, as has been discussed elsewhere“*, it follows that 


KZ 


Now since K lies on the acetic acid axis Wx is zero, so from F 


Equations (2) and (3) 


. 
&x + &z Wz 
By a similar argument it can be shown that for points Y on the 


furfural-rich side 


aad (5) 


0.25 


Ww 
54 66 78 90 


F Ww 
&z &y OR 
——— and ———— : . 
&xk + &z &i + gy Figure 3—Working plot for tie-line determination. 


It follows that experimental determination of 
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TABLE 1 
T1E-LinE DaTA AT 25° (IN PER CENT BY WEIGHT) 





Tie-line 
error 


Water rich phase 


| 
acetic (See 
Water meinaned acid text) 


Furfural rich phase 


acetic 
Water | furfural acid 





| 
| 


87.2 

82.1 | 
76.4 | 
69.2 | 
59.2 | 








is sufficient to fix a tie-line, sincz these quantities are identical to 
respectively Wy/W, and Fy/Fy. A knowledge of the values 
of these latter functions leads directly to the corresponding 
values of W, and Fy for a pair of conjugate solutions, using 
Figure 3. 

As an illustration of the method the following example is 
given for a pair of conjugate phases at equilibrium: 


gz = 18.89 g. 


Wy _ 0.534 
We 


gx + gy 


Thus from Figure 3, Wz is 69.2%, and from Figure 2 the corres- 
ponding values of the haleiad and acid concentrations at this 
point on the binodal curve are respectively 18.1 and 12.7 %. 
Again using Figures 3 and 2, Fy is 71.9%, and the corresponding 
values of the water and acid concentrations at this point on the 
binodal curve are respectively 17.1 and 11.0%. 


Five tie-lines were determined in this manner. The results 
are given in Table 1. In order to determine the accuracy of the 
method the compositions of the synthetic samples in the two- 
phase region used in determining conjugate phase compositions 
were recorded. For absolute accuracy each tie-line should pass 
precisely through the point representing the composition of the 
synthetic sample used in determining the tie-line. The amount 
whereby the tie-line missed each such point at its closest approach 
is also given in Table 1 in terms of % by weight of acetic acid. 
This is a convenient method of describing the quality of the 
data in the present system since the tie-lines are essentially 
parallel to the water-furfural axis. The average amount by 
which a tie-line misses passing through the synthetic sample 
point is slightly less than 0.1% acetic acid. The agreement is 
sufficient to establish that the method is practical. In practice a 
knowledge of the composition of the synthetic samples would 
not be necessary in order to determine the tie-lines by this 
method, but the slight amount of extra labor involved appears 
to be worthwhile when one considers the means it provides for 
checking one’s accuracy. 


It should be noted that for best results points K and L must 
be selected so that in terms of the consolute liquid they contain 
somewhat less than the maximum concentration of that compon- 
ent on the binodal curve. Otherwise for points near the top of 
the binodal curve the addition of neither K nor L will produce 
a locus of points passing through the two-phase region. The 
latter, of course, is essential to the method described here, 
Another point of interest is that for some systems where the 
plait point is displaced considerably away from the maximum 
of the binodal curve, it may be necessary to titrate both of the 
conjugate phases with either K or L for tie-lines near the plait 
point. Such behavior may be anticipated by determining first 
those tie-lines near the base of the two-phase region and noting 
the trend of the tie-lines as the maximum of the binodal curve 


is approached. 


It will be apparent from Figure 3 that the sensitivity of the 
method is not uniform over the entire range of the binodal curve. 
As Fy/Fy or Wy/W, approaches unity the sensitivity decreases. 
Indeed, for any two working solutions such as K and L there 
will be a point of tangency to the binodal curve for each of the 
lines KNZ and YML. At such a point the method of this 
paper is not applicable. The region of the binodal curve subject 
to this difficulty is small, however, judging from the results 
presented here. 
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Nomenclature 


K, L, M, N, Y, Z = various points on the phase diagram (Fig- 
ure 2), either on the binodal curve, the 
furfural-acetic acid axis or the water-acetic 
acid axis; specific points on the binodal 
curve are indicated by the subscripts 
4 2 OF G. 

the weight per cents of water at points N and Z 

on the binodal curve. 

the weight per cents of furfural at points M and 

Y on the binodal curve. 

2x, Zr, Zy, Zz = the weight of samples of solutions wee 

respectively to compositions K, LZ, Y or 


Wy and Wz 
Fy and Fy 
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Machine Computations for 
Transient Stagewise Processes’ 


D. YESBERG? and A. I. JOHNSON® 


A review of recent work on transient stagewise 
mass transfer processes is presented, with particular 
emphasis on numerical and analog solutions for the 
partial differential difference equations involved. The 
numerical or analog solutions are relatively simple, 
whereas the analytical solution even of the linearized 
equation is quite involved. 

An example of gas absorption is considered in 
which the equations are linearized, and solved 
numerically as well as by the use of a resistance 
network analog. These particular methods of solu- 
tion had not previously been used for this problem, 
and the present work is compared with previous 
solutions. The novel feature of the method of solu- 
tion is the representation of the time derivative by 
a backward finite difference to produce a set of 
simultaneous algebraic equations, which may be 
solved most easily by matrix inversion programmed 
on a digital computer. 


r has been pointed out by many authors that a great deal of 
research remains to be done before completely automatic 
chemical plants become a reality“. One of the tools available 
for the study of problems in the instrumentation of automatic 
plants is the ‘electronic analog computer. Using this instrument, 
process equipment may be simulated electrically and the con- 
trol characteristics of such process equipment investigated in 
the laboratory. Besides being useful for the study of control 
problems, an electrical analog of such mass transfer equipment 
as distillation, extraction or absorption towers is useful for pre- 
dicting the time of approach to steady state such as occurs at 
the start-up of the equipment. 

This paper discusses a resistance network analog for a gas 
absorption tower, and its use for predicting the rate of approach 
to steady state. The analog was developed by Liebmann® for 
solving the heat conduction « or diffusion equation. 


Transient Stagewise Processes 


In 1940 several authors discussed transient mass transfer in 
batch fractionation “ 4 > 6 7, However it was not until 
Marshall and Pigford published their book in 1947‘ that the 
concept of transient behavior in stagewise processes was thor- 
oughly developed. Lapidus and Amundson extended the 
special case considered by Marshall and Pigford to a complete 


1Manuscript received October 21, 1959; Accepted January 26, 1960. 
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solution of the transient and unsteady state for a general linear 
equilibrium relation such as is found in absorption and extrac- 
tion problems, and for constant feed rates. The use of a linear 
equilibrium relation restricted the application of the derived 
formulae to limited ranges of concentration changes, but the 
authors felt that qualitative conclusions for the non-linear cases 
could be drawn from the results obtained. Their equation, 
obtained by a mass balance on the m th stage was 


 O¥m yy 4 2Xm = L X,-1(t) + X(t) + 
3 a )=LXa-1 X(t) 


G Yu+i(t) — Y,(t).... 


and by the use of the equilibrium relation 
Y,(t) = a X,(t) + B.. 


results in 
OX n(t) 


(Ha + kh) —— a1 


= LX m-1(t) = 


(L + Gakld) + Cokal..;.......48R 


The conditions at the two ends of the column are described by 
X(t) = f(t) when m = O (upper end) 
Vyii(t) = g(t) when m = M+1 (lower end) 


and the condition of the column at t = 0 is given by 


The Equations (3), (4), (5), (6) make up a complete mathe- 
matical description of the problem, and it is the solution of this 
problem which must be found. 

In 1951, Pigford, Tepe and Garrahan“® solved the equa- 
tions similar to as applied to batch rectification using a mech- 
anical analog computer at the University of Pennsylvania, and 
did not need to limit themselves to a linear vapor- liquid equili- 
brium relation. 

In 1953, Rose et al.“” published a detailed summary of 
their work over the preceding 13 years on the theory of batch 
rectification. Also in 1953, Acrivos and Amundson“ published 
solutions for the problem Lapidus and Amundson had first 
considered in 1950. This time“® a Reeves Electronic Analog 
Computer (REAC) was used to find solutions for both linear 
and non-linear vapor-liquid equilibrium relations. Acrivos and 
Amundson “* considered yet another method of solution using 
matrix methods. 

In 1955, Rose and Williams “» used a REAC to solve the 
problem of transients in continuous distillation, for constant 
relative volatility. They also introduced the dynamic character- 
istics of the stagewise system into the solution, and showed how 
these solutions could determine the most satisfactory method 
of instrumental control. This investigation was continued in 
a later paper “®, 
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In 1956, Jackson and Pigford “” considered the rate of 
approach to steady state in a distillation column. These authors 
replaced the finite difference terms of Equation (1) by the 
approximate differential quotient equivalents, and pointed out 
that the resulting partial differential equation was of the heat 
conduction type. An empirical correlation for the time of 
approach to steady state was derived using the results of com- 
puter calcuations. 

In 1957, Mickley, Sherwood and Reed treated several 
transient stagewise processes in their book “*. A particularly 
important contribution was made by Rosenbrock (9) who over- 
came the disadvantage of the finite difference method of solu- 
tion of Rose et al. “” by using a method similar to that developed 
by Crank and Nicolson @ to solve the heat conduction equa- 
tion. Rosenbrock used a digital computer to solve the equations 
and analyzed his reasons for using this equipment. He also listed 
clearly the idealizations and assumptions that he included in 
the computer program. 

Rose, Sweeney and Schrodt ° have presented an extension 
of the work of Rose et al “” to include multicomponent con- 
tinuous distillation calculations. However, as mentioned above, 
the feasibility of this method for predicting the rate of approach 
to steady state depends on the correct selection of the size of 
the time interval. 

There have been a number of articles published concerning 
the transient behavior of packed columns, but these will not 
be considered. 


Machine Solutions 


The first machine solutions of the partial differential dif- 
ference equation 


OC’ — C'a-t — 2C'n + Cott 
TT. (Ax’)? 


(7) 


where (7) is written in terms of reduced length and time vari- 
ables x’ and T was solved by Hartree and co-workers “) using 
a mechanical differential analyzer. They replaced the second 
order space derivative by the approximate finite difference ex- 
pressions and then proceeded to solve each of the set of ordinary 
differential equations so formed. This method is still used with 
electronic differential analyzers, and a discussion is presented 


by Karplus (23, ch. 12) 


Acrivos and Amundson “*) used an electronic differential 
analyzer to solve ® using circuits very similar in principle to 
those used for the somewhat simpler diffusion equation 
approximation “), 

Another method of machine solution was suggested by 
Beuken in 1937 °, where a direct electrical analogy of the 
diffusion equation may be obtained by considering the flow of 
current along a lumped cable. The electrical equation is derived 
in many texts e.g. °», and for zero inductance, and zero leakage 
to earth from the lumped cable, the equation is 


OV 1 [Vas — 2Vu + Vents (3 
ot RC (Ax)? ria snn Lee 





where R is the resistance, and C is the capacitance per unit 





Figure 1—Lumped cable. 


length of cable. A lumped resistance-capacitance network or 
cable is shown in Figure 1 for comparison with the Liebmann 
circuit to be discussed in this paper. 

Paschkis °® describes the construction of such an electrical 
analog. By a suitable choice of scale factors the electrical model 
may be operated at convenient voltages and time intervals and 
may be built with feasible magnitudes of resistance and cap- 
acitance. The possibility of being able to change the time scale 
is of paramount importance for the practicability of the method, 
A process taking days or hours can be condensed to a few min- 
utes, and vice versa. Paschkis and Hlinka ©’) have extended the 
versatility of this direct analog by the use of unity gain cathode 
followers (3. #: 256), 

Liebmann ® has devised a resistance network analog in 
which the time derivative of “ is also approximated by finite 
differences.  Liebmann’s method of using backward dif- 
ferences °*) satisfactorily overcomes the criticism levelled at 
forward difference representation as used by Rose et al “! ®), 
Liebmann’s analog is discussed in great detail by Karplus 


Time Derivative Approximation 


For numerical solution of the partial differential difference 
Equation (3) the time derivative may be replaced by three dif- 
ferent finite difference expressions. 


Pe 
car 
al 


(X,, is the concentration on the m th tray at time t = nAt.) 
Expressions (9) and (10) are forward and backward differences 
respectively. They have a somewhat larger truncation error for 
a given finite time interval than the Crank and Nicolson expres- 
sion “) @ as mentioned by Crank @ > 18? 


In 1928, Courant, Friedrichs and Lewy “° called attention 
to the fact that in the numerical solution of equations such as 
(7) the choice of the size of the space interval Ax placed a 
restriction on the size of the time interval At when the forward 
difference expression (9) was substituted for the time derivative. 
This restriction does not apply to expressions (10) or (11). 
Von Newman @° *) has developed a method for determining 
whether this restriction applied to other formulae, while more 
elementary discussion of expressions (9) and (11) has been 
published by Rosenbrock “*, and Liebmann °) has discussed 
expression (10). 


Ames a 


Rose et al “':*) have used expression (9) extensively. 
Values of X,,,.+1 are determined successively by repeated appli- 
cation of the equation 


H h 
( ae (Xm a Knn) 


(L + Ga)Xas + Gakatre...........- 08 


to each stage, where the values of X,, at time nAt are known. 
However, as Rose has repeatedly pointed out, the size of the 
time interval is critical for stability of the solution. 


a ED Kiwnte — 


Rosenbrock “® has overcome the instability problem by 
using expression (11) to replace the finite difference on the left 
hand side of (12), while the present work uses expression (10), 
the backward difference, to give a set of simultaneous algebraic 
equations in X,,, (m = 1...M) as unknown, one equation 
for each stage. 


Ha + h 
—— (Xn,.n — Xma-i) = LXu-t9n — 
Al ( n—1) 1, 


(L + Ga)X mn + GaXm+in 
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Figure 2—Resistance network. 


These may be rearranged to 
Xw-in —(1L FO AD) Xan + OX miiin + OXmn-1 = O....(14) 


G H h 
where 6 = ~ and @ = - ie Bd: * (15) 


Equation (14) may be solved using a resistance network, each 
solution giving the values of X.,, for all stages at successive time 
intervals. The analog solutions may then be compared with 
numerical solutions. 


The Resistance Network Analog of Liebmann 
Consider a chain of M equal resistors R,, as in Figure 2. 
Let the beginning of this chain be Po, the first junction point 
be Pi, etc., and the end point be Py. The resistors will be called 
“network resistors”. To each junction point P;, ... Py, ...Py-1 
a resistor of value Ro,, is connected, the free ends of these resis- 
tors, called “series resistors”, being the points Fe eos 
...Py-1. Let the value of Ro,, be related to Ry, by the formula 


feo 
Rint 


The chain of resistors R,, may then be considered a model of a 
stagewise operation. The junction points correspond to stages 
m =0...m = Mim the real equipment, and Pp and Py cor- 
respond to the inlet streams. Let Po and Py be connected to 
voltage sources Vo and Vy corresponding to the prescribed 
inlet liquid concentrations Xp and Xy, and let voltages Vi,o, 

. Vn,o corresponding to the prescribed initial concentrations 
X10, ---Xm.o be applied to the terminals (“feeding points’’) 
Pi, ...P.,... of the series resistances Ro,,. Then the voltages 
Via, ».- Vina, ... Will appear at the network junctions Pi, . . . 
P,, ... and these voltages will correspond to the concentrations 
Xi, ...Xm,1 ... in the medium represented by the model at 
the time t = At. 

If the network voltages V,,,, are measured and the voltages 
at the feeding points P), made equal to V,,,:, then the voltages 
V2 will appear at the network junctions corresponding to 
Xn,2 at the time t = 2At?, and so on. 

To prove this statement, consider the section of the network 
between the points P,,-, and P,,+1 as shown in Figure 3. 

According to Kirchoff’s Law 


3 


ie he, Om oo x nokia (17) 
r= 1 


where /,, I:, 3 are the currents flowing to the junction points 
P,,. Thus 


V (eee eer). 4 
m—\,n Rms Re. mn 


Re V, + Xn V, =0O (18) 
Rint sai Rom er ey cae” ah : 


which is formally identical with (14) provided (16) holds. 
Hence the voltage relations in the described network model 
represent the difference Equation (14) within close limits set 


by the accuracy of the voltage measurements and the tolerances 
of the resistance values. 
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Figure 3—Liebmann analog element. 
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Figure 5—Comparison of analog and numerical solutions 
and illustration of the h? extrapolation. 
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The solution of the problem by the described analog model 
is subject to truncation errors of order At, but it is stable for 
any value of At, unlike solutions utilizing expression (9). 


Apparatus 


An apparatus similar to that of Liebmann ® providing 11 
network points was constructed, with provision made for use 
as two sets of six network points. It was assembled in two relay 
racks 44 inches high as shown in the photograph Figure 4. 
The relay rack on the left contained the resistance network, 
while the relay rack on the right contained the potentiometers, 
series resistors, and a set of resistors for allowing leakage to 
earth from the network points—-“the leakage resistors’. Further 
details may be found in Yesberg “. The use of the equipment 
may best be evaluated by considering a suitable numerical 
example. 


Example of Gas Absorption 


As a simple numerical example, Lapidus and Amundson (9) 
considered a six-plate column with solute component A in air 
and with water as the absorbent. The fat gas composition jumped 
from 0.2 lb. A per lb. of air to 0.3 Ib. A per ‘Ib. of air. The problem 
is to determine how the fat solvent and lean gas compositions 
vary with time, and how long it takes to reach within 1% of 
steady state. Table | gives the required numerical data. 


TABLE 1 
_ Gas ABSORPTION EXAMPLE 





Numerical Value 


0.72 
0.0 
66.7 


Symbol 





_ 
= 
oo 


75 
1 


L 
Ga 


on 


M 

X, for all ¢ 
Yur, t = O 
Yau, i > O 
d 


cee. SY 
MmwBwnNnodc Go oo 


~ 
~ 


Lapidus and Amundson “ first found the steady state 
values of A,, corresponding to Yy+1 = 0.2 by use of the steady 
state form of their solution to Equation (3). Table 2 gives these 
values of A,,. 

TABLE 2 
INITIAL STEA ADY STATE SoL UTION 


With these steady state values known, the variation with 
time of X the fat oil composition and Y; the lean gas composi- 
tion can be found. It should be pointed out, however, that some 
transient values as presented in Lapidus and Amundson’s 
Table 1 “ were in error; but these errors were corrected in a 
later paper “®. The correct values are presented in Table 3. 
The steady state fat oil and lean gas compositions can be shown 
to be 0.3821 and 0.0662 respectively. 


Analog Computation 


Four computations were performed on the analog to deter- 
mine values of X for each tray ten minutes after a step change in 
the concentration of the entering fat gas stream had occurred. 
The time of ten minutes was chosen because it was known that 
the truncation errors would be significant at this time. The 
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numerical coefficients for the sets of simultaneous equations, 
which determine also the ratio of values of adjacent resistors ip 


the analog, are recorded in Table 4. 


TABLE 3 
SOLUTION OF Gas ABSORPTION Ex AMPLE 


Time 











TABLE 4 
Nu ME RICAL CoEF FFICIENTS 





Ga 
Intervals | Minutes | L 


477 
.177 
RTT 
S77 





Details of each of the first two computations are presented in 
full, the results of all four are plotted on Figure 5. 


Notice that (14) and (15) are satisfied by the value of the re- 
sistances. The values for the initial steady state depend only 
on the values of the resistances R,, and the voltages applied at 
the ends of the network R,,. The difference between the expected 
values and the measured values is an indication of the accuracy 
of the component resistors and efficiency of the switches and 
potentiometers. Since the voltage across the potentiometers was 
approximately one volt, the voltage errors are of the order one 
to three millivolts. There are two concentrations of interest, 
the liquid composition on the bottom of the tray, and the gas 
composition leaving the top tray. For the purpose of comparing 
the numerical values, however, only the values of liquid com- 
position on the bottom tray are presented. 


Numerical Solution on a Digital Computer 


To find the expected values of the voltages on the analog, 
numerical solutions of the sets of simultaneous ee (14) 
were performed on an IBM 650 digital computer. A slight 
variation of the Fox and Thomas modification of the Choleski 
method of matrix inversion “: #1) was programmed using 


the S.O.A.P. II procedure “. 


(a) Steady State Solution 


To find the initial steady state solution At was set equal to 
infinity i.e. @ = O in (14), and the set of six simultaneous 
equations of the form 


-( ) 3 x Ga 
Xe —(14+—)X+— 


was solved for the six unknown concentrations differing by at 
most one in the fourth decimal place from the values in Table 2. 
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TABLE 5 
ONE TIME INTERVAL* 





R, (ohms) | 569 





| | 
R, (ohms) 3610 








Stage No. 





Initial 
Steady State 





Expected 
Measured 
Difference 














TABLE 6 
Two TIME INTERVALS* 
= = — = — = = — = T = = =T — 
R, (ohms) 1340 














R, (ohms) 4248 








Stage No. 








Initial 
Steady State | a x ; oe ' .455 





Expected ; 4 «Sor 
Measured iG. , é ‘ : . 554 
Difference ‘ 0 | — .003 




















-611 
— .001 


Measured x 4 390 | .501 
Difference | .001 —.001 


| “7 
Expected = | .391 502 | 612 
| 


*Values of the voltages are greater than the values of the corresponding concentrations by a factor of two. 


Since there is no truncation error involved in Equation (19), it TABLE 7 
Is an exact expression. APPROACH TO STEADY STATE, ¢ = 10 MINUTES 
The same method may be used for the final steady state. : ee ] 





| | | 

tray tray | tray tray | tray tray 
| 2 3 4 5 6 
| | 

1346 | .1911 | .2444 | .2974 

step change in inlet gas concentration occurred, until one hour .07284 | .1368 | .1953 | -2511 | .3063 

after the change, when the concentrations were easily within “ -07354 | .1385 | .1984 | .2556 | .3117 

one per cent of steady state. Time intervals from 0.25 minutes ‘oo tr ‘a oa ‘ia 

to 60 minutes were used. ae | “4408 | | oe: * 

Table 7 presents the results at ¢ = 10 minutes on each trav. 
These values are used to draw the continuous lines in Figure 5, 
on example of f?-extrapolation (35). An analytical value was 
available for tray 6 from Acrivos and Amundson’s work, and 
it is apparent that h?-extrapolation is effective in this case, using 
either the analog results or the exact numerical solution. 


(b) Transient Solution At mins. 


The computer program was arranged so that the concentra- 
tions on each tray were calculated from time zero, when the 07186 





| 
| 





APPROACH TO STEADY STATE: BOTTOM TRAY 


| | 
, 3 : a Time Lapidus & t = 0.25 min. 
Table 8 is an alternative presentation of the digital computer unten. | Amundson 


results, and these results show that for a finite time interval 
of 0.25 minutes, the numerical computation is within one per | _ 
cent of the analytical values for tray 6. Figure 6 is a plot of ‘oa 
these data together with the results of the digital computer com- | 3443 

| 





putation for a finite time interval of five minutes. 3562 
. 3697 
.3754 
.3785 


Discussion 


A severe criticism of the Liebmann analog as presently con- ti 
structed was the time required for the solution of a problem. . 3821 


While it was very much faster than any manual numerical method ln Ses i eee 
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Final steady state 


Lbs. solute per Ib inert liquid 


—— Exact numerical solutions 
© Analytical solution 
( Acrivos & Amundson ) 


20 30 
Time in minutes 


Figure 6—Rate of approach to steady state of liquid leaving 
the bottom tray. 


of a similar accuracy, it was no faster than programming a 
digital computer to carry out the computation. The British 
Admiralty has had a Liebmann analog constructed in which 
potentiometers are balanced by servo mechanisms “®, This 
modification would speed the solution enormously. 


A second criticism is the fact that this analog cannot be 
used in the investigation of control dynamics. If an analog 
of the Paschkis type had been used, the truncation errors 
would vanish for the solution of the problem on gas absorption 
in a stagewise process and the solution would proceed auto- 
matically, thereby being extremely useful for investigating con- 
trol dynamics. However, the cost of the Paschkis type of equip- 
ment would be far greater and the experimental troubles of 
greater complexity, particularly the introduction of non- 
linearities as the solution proceeds. Non-linearities should be 
easily introduced on the Liebmann analog, for example by vary- 
ing a in Equations (13) or (14) with concentration; however, 
the authors have not investigated this modification extensively. 

It is suggested by this work that the use of forward finite 
differences to approximately represent a time derivative might 
be discontinued in favor of either backward differences following 
Liebmann’s work, or by using the method of Crank and Nicolson. 
Perhaps the backw ard differences are simplest when a digital 
computer is available. 


‘ 


Nomenclature 


composition of liquid on the n th plate at time ¢ = O; 
lb. solute 1 Ib. inert absorbent 

= capacitance; farads 

= concentration of diffusing material 

= boundary condition of liquid flowing to the top plate; 
a function of time 

boundary condition of gas flowing to the bottom plate; 
a function of time 

inert gas rate; lb./hr. (solute free) 

vapor hold-up; Ib. of inert gas 

liquid hold-up; Ib. of inert liquid 

electric current; amperes 

absorbent rate; Ib./hr. (solvent free) 

stage number 

total number of stages 

time interval number 

resistance; ohms/unit length 

series resistor; ohms 

network resistor; ohms 

dimensionless time 

time; hours or minutes 

electric potential; volts 

concentration of liquid leaving the m th plate; lb. solute/ 

Ib. inert absorbent 

space variable; length 

dimensionless space variable 

composition of vapor leaving the m th plate; lb. solute 
inert gas 


Ds BeOS a 
3 


3 


STN 
2 


| 





vapor liquid equilibrium constant 
vapor liquid equilibrium constant 
finite difference operator 

Ga 

L 

Ha —h 


~ LAt 
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Preparation of Aluminum Octanoate and 


Aluminum Isooctanoate’ 


J. A. WHEAT?, D. A. MacLEOD?, R. G. DUNLOP?, 


and H. 


Two aluminum soaps have been prepared as 
hydrocarbon gelling agents. The aluminum soap of 
2-ethylhexanoic acid (Octal) was successfully pre- 
pared by a pilot plant continuous process which 
utilized a flash mixer. Optimum process conditions 
and the effect of important process variables were 
determined. There was a seasonal cycle of short 
gelling times and low gel consistencies in the summer 
months and high results in the winter. The cause of 
the cycle is not known. 


Aluminum isooctanoate (Isooctal) was prepared 
continuously and batchwise on both a laboratory and 
pilot plant scale. The continuously prepared products 
were not satisfactory. Good Isooctals were prepared 
batchwise and pilot plant results were the same as 
those obtained in the laboratory. Gel consistencies 
showed a seasonal cycle but thickening rates did not. 


an soaps have been widely used to thicken 
hydrocarbon solvents in the manufacture of lubri- 
cating greases, paints, inks, and cosmetics. During World 
War Il their use was extended to the thickening of gaso- 
line for incendiary bombs and flame throwers. This paper 
describes some of the work done in these laboratories on 
methods of preparation of aluminum soaps for this latter 
purpose. 

Aluminum soaps for thickening hydrocarbon solvents 
are prepared by combining in aqueous solution the fatty 
acid, alkali, and an aluminum salt. They consist mainly 
of aluminum hydroxy disoap but also contain small quan- 
tities of aluminum dihy droxy monosoap, free fatty acid, 
and aluminum hydroxide “2%. 4). Aluminum trisoaps 
have been prepared (%: 4) under anhydrous conditions but 
are not formed in aqueous solutions nor can they exist 
in the presence of water. The overall equation of the 
formation of an aluminum disoap is 

Al,(SO,), + 4RCOONa + 2NaOH 
— 2Al(OH)(OOCR), + 3Na,SO, 
The reaction mechanisms have been studied (5) but will 
not be discussed here. The theoretical alkali requirement 
for disoap formation is seen to be 50% in excess of that 
required to convert the fatty acid to the sodium soap. 
This excess alkali is usually sodium hydroxide but sodium 
carbonate has been used in the preparation of aluminum 
stearate (6), Aluminum sulphate is the most commonly 
enjoyed aluminum salt but the chloride can also be used 7). 
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SHEFFER? 


In the simplest batch production method, a solution 
of aluminum sulphate (alum) is added to a sodium soap 
solution containing the excess sodium hydroxide, or the 
alkali-sodium soap solution can be added to the alum. 
Another variation is to add the alum and excess alkali 
simultaneously to a sodium soap solution. In general, there 
are a number of additional ways in which three reactants 
can be combined either simultaneously or in sequence but 
in this particular case most of them cannot be used because 
mixing the alum and excess alkali would produce insoluble 
aluminum hydroxide. 


To combine liquids continuously, equipment may 
range from flash mixers such as a spinning disc or a 
centrifugal pump to stirred tanks either singly or in series 
(8). Stirred tanks may follow flash mixers to give increased 
retention time or to serve as surge tanks; they may be 
operated continuously or batchwise. Woods and Tay lor (+ 
developed a multiple stirred tank system for aluminum 
soaps and a British patent ‘®) describes a similar system 
but with the alum added in two stages and the excess 
alkali added to a third vessel. In these laboratories, con- 
tinuous methods based only on flash mixing devices were 
investigated. The essential difference between any con- 
tinuous process and the simple batch process is that in the 
former all the product is formed at a single pH level 
while in the latter the pH changes from about 12.5 to 
3.5-4.5. A continuous process based on a flash mixer differs 
from stirred tank continuous processes as well as batch 
processes in that the product is quickly withdrawn from 
the reaction zone and not a'lowed to contact unreacted 
chemicals. More reproductible products would therefore 
be expected from the flash reactor system. 

The product is separated from the mother liquor by 
filtering or centrifuging and, after washing with water, 
is dried and ground. Most soaps are easily filtered or 
centrifuged but washing is quite slow because sulphate 
is strongly held by the soap. Prolonged w ashing is detri- 
mental because hy drolysis is likely to occur. It is impor- 
tant that the wet material be comminuted before drying. 
It is generally believed that dryi ing is a critical step and 
high speed, low temperature dryi ing as in a flash dryer 
is favored but rotary ‘!®) and tray dryers can be used. 


Gasoline is thickened by mixing with it the desired 
amount of aluminum soap; with some soaps heating is 
necessary and frequently a peptizer is added. The peptizer 
increases the rate of gel formation but also increases the 
rate of deterioration of the gel. Some peptizing action is 
obtained from residual water and free fatty acid but this 
is usually not sufficient. The more common peptizers 
include alcohols, Cellosolves, phenol, cresols, xylenols, 
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Figure 1—Continuous processes used to prepare aluminum 
soaps on a pilot plant scale. 


and fatty acids ‘471!!), The polymeric structure and 
properties of the gels are described in a number of papers, 
for example (3. 4,7, 11, 12), 

An ideal aluminum soap for thickening gasoline 
would have the following properties: It would readily 
dissolve in gasoline and rapidly form a gel over a wide 
range of temperatures with a minimum of stirring; A low 
concentration would give a high viscosity; The soap and 
gels prepared from it would be stable for long periods 
of time under a wide range of temperature and humidity 
conditions, The soaps and gels would be reproducible; 
and the soap would have a high bulk density. 


It has been well demonstrated (1,13) that aluminum 
soap-gasoline gels possess properties that make them good 
incendiary fuels but this aspect was not considered in the 
work being reported here. All soaps were compared on 
the basis of laboratory tests only. 


The first aluminum soaps used extensively to thicken 
gasoline for fire bombs and flamethrowers were aluminum 
stearate (1!) in the United Kingdom and Napalm (1°, 18, 14) 
in the United States. Napalm was made by coprecipitating 
a mixture of naphthenic acids, oleic acid, and the acids 
from coconut oil. The principal constituent of coconut 
oil is lauric acid and this led to the study of aluminum 
laurate (4%), A number of other fatty acids such as mixed 
isooctanoic acids (15) and 3,5,5-trimethylhexanoic acid (16) 
have been investigated. In these laboratories, experimental 
work. was concerned largely with four acids, 2-ethylhexa- 
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noic acid, mixed isooctanoic acids prepared by an oxo (1) 
process, stearic acid, and lauric acid. Data for the first 
two of these are reported in this paper. 


EXPERIMENTAL 


Soaps were prepared batchwise and continuously on 
both a laboratory and pilot plant scale. Most of the work, 
however, was done on a pilot plant scale. The laboratory 
processes were introduced because supplies of isooctanoic 
acid were limited. 

Figure 1 illustrates the four continuous or semi- 
continuous processing methods used. In Method A, sodium 
soap plus excess sodium hydroxide was combined with 
alum solution and the resulting slurry held for a time in 
a tank before centrifuging. In Method B, the slurry was 
sent directly to the centrifuge. In the other two methods, 
sodium soap and alum solution were combined in a flash 
mixer; in Method C excess alkali was added simultaneously 
to the slurry tank while in Method D the alkali was added 
batchwise after stopping the flow of the other two 
solutions. 


Pilot Plant Equipment 


For pilot plant continuous runs 10 lb. of fatty acid and 
the desired amount of sodium hydroxide were dissolved 
in water in one of the solution tanks shown in Figure 2 
to give 22 gal. of solution. Alum solution was prepared 
in the other tank so that an equal volume gave 5-10% 
excess alum. These solutions were pumped at equal rates 
of about 0.75 g.p.m. by the two halves of a duplex 
reciprocating pump. The pulsations in the discharge were 
smoothed out by a system of orifices and surge tanks. An 
orifice diameter of about 1/16 in. was used and the surge 
tanks were 3 in. in diameter by 15 in. long. The back 
pressure built up to 20-30 p.s.i. gauge and fluctuated only 
2-3 p.s.i. over the suction-discharge cycle of the pumps 
so the liquid flow from each orifice was essentially smooth 
and uniform. When excess alkali was added separately it 
was pumped with both halves of a second duplex recipro- 
cating pump with a total capacity of 4.9 g.p.h. Pulsations 
in the flow of this third solution were not detrimental. 
The alum and sodium soap solutions were led by means 
of concentric pipes to the inside of either a colloid mill 
or a centrifugal pump. The colloid mill, which was used 
first, had a 3 in. diameter rotor driven at 7300 r.p.m. 
Since it was usually operated at maximum clearance of 
0.05 in. any spinning disc would probably have worked 
as well. The centrifugal pump had an open 6 in. impeller 
and a speed of 3450 r.p.m. Both units had 3 h.p. motors. 
The flow of solutions was stopped after 24 min. and the 
remaining processing steps carried out batchwise. 


Slurry was transferred to a centrifuge using a centri- 
fugal pump and the product collected in a 26 in. diameter 
basket operated at 855 r.p.m. It was washed at 1-2 g.p.m 
for about 20 min. and then spun for 5-10 min. at 1710 
r.p.m. It was then removed with an unloading blade, 
spread 0.1-0.5 in. deep on screen bottom trays 18 x 28 x 2 
in. deep, and placed in a tray dryer. The dry bulb 
temperature of the dryer was automatically controlled by 
regulating the flow of steam to a heater and the wet bulb 
temperature by regulating a steam spray. The air velocity 
was about 250 ft./min. A constant drying time of 21 hr. 
was used and, except for special tests, Octal was dried at 
60°C. and Isooctal at 78°C. A higher temperature was 
needed for drying Isooctal not to increase the rate of 
drying but because Isooctal had a higher equilibrium 
moisture content. Dried soap was removed from the dryer 
and ground in a 6. in. hammer mill. 


The same equipment was used for batch precipitation 
by preparing the sodium soap solution in the slurry tank 
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Figure 2— Diagram of pilot plant equipment used to 
prepare aluminum soaps. 


and adding a relatively concentrated alum solution via 
the excess alkali system over a period of about 18 min. 

All equipment and piping except the copper wash 
water tank were 304 or 316 stainless steel. 


Laboratory Apparatus 

Batch and continuous laboratory runs were made with 
50 g. of fatty acid in 500, 1000, or 2000 ml. of solution 
and with a reaction time of 18-20 min. In batch operation 
the reactants were mixed in a beaker by a 2 in. propeller 
with a variable speed drive. The rate of flow of alum 
solution was set by a constriction which controlled the 
admission of air into a separatory funnel. The product 
was recovered in a variable speed, 5 in. centrifuge and 
dried on 8 x 8 x 4 in. trays in the pilot plant dryer or 
a small electrically heated dryer. 

In continuous operation only Method B was used with 
the reactants being mixed in a small centrifugal pump. 
At first, flow rates of the solutions were regulated with 
needle valves and measured with calibrated rotameters but 
a peristaltic action pump (Sigmamotor, Inc. Model T-63) 
installed later was better. 


Materials 
Technical grade 2-ethylhexanoic acid and mixed iso- 
octanoic acids were purchased in 400 lb. drums from 
Carbide Chemicals Company, Montreal but were produced 
in the United States. The isooctanoic acid was the oxi- 
dized product of an oxo process ‘!7) and the following 
data were supplied. 
Sp. gr. 20/20°C. 0.9171 
Distillation at 760 mm. 
Initial b.p. 229.0°C. 
50% 231.7°C. 
End point 239.2°C. 
Purity 96.9% 
Water 0.04% 
Ester 0.08% 
GCdor 70° 


Flake caustic soda was purchased in 380-400 Ib. drums 
and commercial grade, iron free aluminum sulphate 


(Al,(SO,), . 18H,O) was purchased in 100 Ib. bags. 


Test Methods 

Octal gels were prepared by dispersing 1.5% Octal in 
gasoline with 0.5% distilled 2-ethylhexanoic acid as pep- 
tizer. Isooctal was tested at a concentration of 2% with 
no peptizer. The time for the vortex to disappear when 
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mixing the ingredients under standard conditions was 
reported as stir time. Gel consistency was measured with 
a Gardner mobilometer and the results were reported as 
the weight required to give a constant rate of shear. Gels 
were tested first after two days and then each succeeding 
week from the date of preparation. These later results 
were reported as a percentage of the 2-day value. Ordinary 
leaded gasoline was used for all gel tests. When replace- 
ment was necessary, samples of gasoline were obtained 
from Ottawa distributors and compared against a standard 
lot of Octal. A drum of the gasoline most closely resemb- 
ling the previous drum was purchased. 


The moisture content of both wet and dry soap was 
determined with an infra-red moisture balance (Central 
Scientific Company, Chicago) with a 125 watt infra-red 
lamp mounted about 1.5 in. from the sample. Octal 
samples were left on the balance until they reached 
constant weight. For wet material, 30-40 min. was required 
while for dried product, 10 min. was sufficient. Isooctal 
samples decomposed if left too long and it was necessary 
to plot loss in weight against time and take the first 
abrupt change in slope as the moisture content. This 
change in slope was usually readily detected. All moistures 
were reported on the wet basis. 

Bulk density was determined by measuring the volume 
of about 35 g. “of dry product in a 250 ml. graduate. The 
cylinder was tapped on the bench top either manually 
or mechanically to obtain a minimum volume. 


RESULTS — OCTAL 


Seasonal Variation 


Throughout the work on this project there were 
several unaccounted for changes in stir time and Gardner 
consistency. Since these appeared to be related to time 
of production, graphs were prepared of weekly average 
results plotted against numbers representing the weeks of 
the year. Figures 3 and 4 show the changes of weekly 
stir time and 2-day Gardner consistency with date of 
preparation for all lots of Octal that could be reasonably 
included in the averages. The same data are presented in 
Figures 5 and 6 after smoothing by means of a seven-point 
moving average. These graphs show a definite seasonal 
cycle that is remarkably reproducible from one year to 
the next. 


Several samples of Octal were tested over a period of 
time and, as shown in Table 1, these results followed the 
same seasonal trend as freshly produced samples. How- 
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Figure 3—Variation of weekly average stir time of Octal 
gels with date of preparation. 
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Figure 5—Lines of Figure 3 after smoothing by plotting 
a moving seven point average. 


ever, when winter and summer soaps were tested at the 
same time, the winter soaps did give higher results than 
the summer soaps. Thus, the test results depended upon 
both time of production and time of testing. Weekly 
average results of 1954-1955 (Figures 3 and 4) were 
compared with weekly averages of daily mean temperature 
and daily mean absolute humidity. As expected, the 
general trend of the curves were in agreement, however 
weekly fluctuations in test results did not correspond with 
Ww eekly fluctuations in temperature or humidity. 


TABLE 1 
EFFECT OF TIME OF TEST FOR Two Groups OF OCTALS 
PREPARED WITH 35° Excess SopIUM HyDROXIDE 


Stir time, 2-day Gardner 

Week of test min. consistency, g. 
1 Nov. 1954* 10 186 
17 Jan. 1955 21 176 
21 Feb. 1955 21 187 
28 Mar. 1955 34 184 
6 Jun. 1955 s 139 
28 Feb. 1955* 68 213 
13 Jun. 1955 14 125 


*Time of Production 
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Figure 4—Variation with date of production of weekly 
average 2-day Gardner consistency of Octal gels. 
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Figure 6—Lines of Figure 4 after smoothing by plotting 
a moving seven-point average. 


The direct cause of the seasonal cyle is not known but 
its presence had to be taken into account in planning 
experiments and interpreting results. In all experiments, 
adequate controls were included and results compared 
with these controls rather than previous data. Factorial 
experiments, which sometimes took from three to six 
weeks to complete, were carried out in randomized blocks 
and nearly always there was a statistically significant block 
effect in accord with the seasonal cy cle. 


Preparation of Solutions 


The effect of excess sodium hydroxide on both summer 
and winter soaps is shown in Table 2. As the amount was 
reduced from the theoretically required 50%, the stir 
time became shorter and initial gel strength higher but 
the gels tended to be less stable. Figures 5 and 6 also 
show that 35% excess alkali gave lower stir times and 
higher gel consistencies than 50%. A level of 35%, which 
was selected as optimum, had the disadvantage that the 
product yield was about 75% of theoretical as compared 
with 80% when 50% excess sodium hydroxide was used. 
All yields were low probably because of incomplete 
reaction or hydrolysis during washing since mechanical 
losses were not large. 
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TABLE 2 


E FFECT ( OF F SOLU TION VARIABLI ES ON N GEL L FORMING PROPE RTIES OF Oct AL 


Excess sodium Number 
hydroxide, 7 % of runs 


30 12 
0.45 35 12 
0.45 40 12 
0. 45 50 12 


Soln. concn., 
lb. acid/gal. 


: Stability, % 
2 day Gardner 
consistency, g. 21 days 


Stir time, 

min. 7 days 
a, 170 47 24 
8 170 54 24 
12 160 66 30 
17 100 82 52 





30 
35 


.23 
45 
91 
45 
on 





45 
45 
45 








Increasing the level of excess alum from 10% to 30% 
had no effect. This means that relative pumping rates of 
the two solutions could vary over a wide range without 
affecting the product. 

As the reactant solutions were concentrated, stir time 
decreased markedly (Table 2) but gel consistencies were 
not affected. Thus, concentrated solutions are preferable. 
This is an economic advantage as well. For most of the 
experimental work however, the intermediate concentra- 
tion of 0.45 Ib. of acid/gal. was used because this concen- 
tration gave the same results as 0.91 lb./gal. but the 
consumption of chemicals was lower. 


There was no difference between reaction temperatures 
of 25 and 40°C. but at 70°C. the stir time was higher 
(Table 2). With 50% excess sodium hydroxide, the tem- 
perature effect was greatly magnified and soaps prepared 
at 70°C. had very long stir times and gave very weak gels. 


Reaction Conditions 


There was no difference between Octals prepared with 
the colloid mill and those prepared with the centrifugal 
pump. A centrifugal pump would be preferred because it 
has a lower first cost. Reducing the flow rates by one- 
half did not improve the product so the higher rate is 
preferable. Data on the effect of mature time and a three 
solution method of preparation are given in Table 3. It 
is apparent that the three solutions method gave an 
inferior product and that 60 minutes mature time was 
neither necessary nor detrimental. 


TABLE 3 
ErrEct OF MATURE TIME AND PRE PARATION METHOD (Sera. ) 








7 | 
| 
Num- | | 2 day 
Reten-| ber | Gardner 
tion of time, | consist- 
time, | runs min. ency, 
min. | g. 


Produc- Stir Stability, 
ton 
method 
7 days 
0 11 192 | 67 | 25 
60 12 194 74 25 


21 days 





0 | 37 | 198 | 69 19 
_ | 40 | 198 72 19 


0 49 | 209 99 62 
0 week aD cam 
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46 200 85 28 
64 206 78 36 
71 177 87 51 
87 128 31 


17 201 4 32 
11 193 24 
10 191 64 26 


3 203 
54 208 
9 169 
9 158 
14 158 


Centrifuging and Drying 

With the standard 10 Ib. of acid per run, the centrifuge 
cake was about 1 in. thick and had a moisture content of 
55%. With only 5 lb. of acid the moisture content was 
65% and when 20 Ib. of acid was used it was 45%. In two 
special runs the centrifuge basket was filled to the maxi- 
mum thickness of 3.5 in. with no detectable decrease in 
filtration or washing rate. The average moisture content 
was about 40%. 


The effect of final moisture content on the gelling 
properties of Octal is shown in the first portion of Table 
4. These data were obtained by dryi ing trays with different 
amounts of Octal for a uniform time of 21 hr. Using an 
average slope, it was calculated that increasing the 
moisture content from 1% to 2% reduced the stir time 
2.7 min. and the 2-day Gardner consistency only 14 g. 
Even at almost 4% moisture, the gel consistency was still 
158 g. 


The second portion of Table 4 shows the effect of 
dryer temperature. The tray loading was also varied in 
order to obtain a more nearly constant final moisture 
content. Under these conditions dryer temperature from 
60 to 80°C. had little effect. At 90°C. gel consistencies 
were lower and at 50°C. stir time and gel stability were 
lower. Some Octals were freeze dried at a maximum tem- 
perature of 24°C. and the average stir time was about 


TABLE 4 


EFFECT OF DRYING CONDITIONS ON GELLING 
_ PRoPE RTIES ( OF F OcTAL 


2 day 
Gardner 
j } Stir consist- | 21 day 
temp., | Someones |Moisture,| time, ency stability, 
i | Ib. /sq. ft. Q% min. g. q% 
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tions were changed, and drying temperature was reduced 
from 78 to 60°C. without improvement. Nor were the 
results improved by the addition of acetic acid to the 
sodium soap solution ‘+ *) or of 2% silical gel to the dried 
product. 



























































> Production by Batch Process 
& As with Octal, 2-day Gardner consistency results 
ti showed a pronounced seasonal cycle but the stir time § ——— 
o cycle, if present at all, was very small. The points in , 
a s, oe No 
= Figure 7 are weekly average stir times and 2-day Gardner Tw 
° consistencies while the lines were drawn through moving Am 
= five-point averages with the weekly averages weighted Ace 
” according to the number of runs available. Only the results om 
for laboratory soaps prepared under standard conditions —— 
£ were used. These standard conditions were as follows: 
7 Sodium soap solution concentration 5% acid _Bo 
= Excess sodium hydroxide 40% with ¢ 
F Alum solution concentration 40% tion tI 
« Excess alum 5% the ra 
o Reaction temperature 40°C. solutic 
40 50 Reaction time 18 min. prope! 
WEEK Mature time 10 min. it 
Figure 7—Effect of date of production on properties of Wash water 40 g./g. acid could 
gels of Isooctal prepared batchwise on a laboratory scale. Drying temperature 76°C. Thus, 
Drying time 21 hr. cake \ 
one-half that of material tray dried at 60°C. The initial Tray loading (laboratory 0.25 Ib./sq. ft. g./lite 
gel consistencies were the same but gels of the freeze : ; (pilot plant) 0.7 Ib. /sq. ft. bariun 
dried product were less stable. Final moisture content 0.6-1.2% only < 
Excellent results were obtained when Octal was flash As seen in Figure 7, stir times were 10-15 min. and aur 
dried in a No. 40 Raymond Imp Mill Flash Dryer at 2-day Gardner consistencies varied from 70 g. in the but 2- 
Mallinckrodt Chemical Works Limited, Montreal. Gel summer to 155 g. in the winter. Gel stabilities were about In 
consistencies were the same as those of tray dried Octals 42% at 7 days and 20% at 28 days. X-45, 
but stir times were reduced to about one-third, although TABLE 5 — 
this may have been largely because the tray dried material EFrFeEct or Process VARIABLES ON GEL ForMING Caen 
was not ground. These data cannot be compared with PROPERTIES OF ISOOCTAL —__ it was 
the results obtained here because the soaps were prepared a (i a ia | 7 were | 
batchwise. Nevertheless, flash drying would be an ideal im ; 2 day except 
method for large scale continuous processing. Excess Reac- | Reac- | Wash | Stir | Gardner alum. 
sodium tion tion | water, | time, | consist- difficu 
Miscellaneous hydroxide, time, g./g/ | min. ency, ow adi 
The average bulk density of Octal was about 0.20 , aun. | acid &. soap s 
g-/ml. and on a Tyler “Rotap” apparatus 98% passed 100 113 and st 
mesh and 13% passed 200 mesh. La time. | 
Summary 131 materi 
’ De " of the 
The optimum or standard conditions for preparing 131 Anal 
Octal by the pilot plant continuous process were as 170 slight! 
follows: 183 jae 
Sodium soap solution 129 
concentration 4.5% acid (0.45 lb./gal.) 170 
Excess sodium hydroxide 35% 172 
Excess alum 10% ees 
Reaction temperature 40°C. 
Pumping rate 0.750 g.p.m. 
Wash water 35 |b./lb. acid ‘ : Ss 
Drying temperature 60°C. ———— a 
Drying time 21 hr. 
Tray loading 0.5 Ib./sq. ft. RINSING WITH CALCIUM AND BARIuM SoLutions (IsooctTaL) Labora 
Final moisture content 0.6-1.2% (Number of runs = 4) Pilot P 
RESULTS — ISOOCTAL va = 
Production by Continuous Process Rinse Stir Gardner 28 day 
All Isooctals prepared continuously were very poor. solution pom consistency, | stability, 
The average stir time for 65 laboratory runs using Method z ; ~ 
B was 32 min. and the average 2-day Gardner consistency 
was 36 g. The averages for 38 pilot plant runs using None 12 
Method A with a mature time of 5 or 10 min. were 100 BAtOHS - 
min, and 75 g. Excess sodium hydroxide was varied from CaCl, : 7 
35 to 50%, reaction temperature was tested over the range BaCl, 9 
from 30 to 45°C., pumping rates and solution concentra- = ca as _ 
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TABLE 7 


EFFECT OF ADDITIVES (ISOOCTAL) 
AS to anne: solution, Number of runs = oe 


Stir Gardner 28 day 
time, consistency, stability, 
min. g. % 


Additive 


None 15 : 16 
Tween 40 13 11 
Ammonyx CO 10 10 
Acacia gum 11 | 12 
Acetic acid 13 | ‘7 


Both stir time and gel consistency increased slightly 
with excess sodium hydroxide, reaction temperature, reac- 
tion time, and amount of wash water (Table 5) but over 
the range from 2.5-10 g. of acid/100 ml., sodium soap 
solution concentration did not affect the gel forming 
properties. 

It was thought that sulphate remaining in the product 
could be precipitated with barium or calcium solutions. 
Thus, following the regular water wash, the centrifuge 
cake was rinsed with 500 ml. of solutions containing 1.5 
g./liter of barium or calcium hydroxide or 15 g./liter of 
barium or calcium chloride. Of the four solutions tried 
only calcium hydroxide looked promising (Table 6). 
Stir time was considerably lower and stability was higher 
but 2-day consistency was not as high. 

In an experiment on the effect of additives (Triton 
X-45, Tween 40, Ammonyx CO, acetic acid, adipic acid, 
gelatine, glue, locust bean gum, guar gum, acacia gum, 
carboxymethyl] cellulose) at a level of 1% of the acid 
it was found that it made no difference whether materials 
were added to the sodium soap solution or to the alum 
except that a number of them were precipitated by the 
alum. On a large scale this precipitate would cause no 
difficulty but on a laboratory scale it made it impossible 
to add the alum solution slowly and uniformly to the 
soap solution. All additives tested increased gel strength 
and stability but nearly all of them also increased stir 
time. In a second experiment, four of the most promising 
materials were tested at levels of 0.4, 1.0, 2.0 and 5.0% 
of the acid. Results for the 1% level are given in Table 7. 
Acetic acid and Ammonyx CO were the best, giving 
slightly lower stir time and considerably higher gel con- 
sistencies. 


Batches prepared on a pilot plant scale gave the same 
results (Table 8) as laboratory batches prepared during 
the same period. 


DISCUSSION 


Octal was successfully prepared by a pilot plant con- 
tinuous process but Isooctal had to be made batchwise. 
However, Isooctal had a bulk density more than 2.5 times 
that of Octal and it could be used without a peptizer. 
Both stir time and Gardner consistency of Octal gels had 
a pronounced seasonal cycle but only the consistency of 
Isooctal gels changed w ith time of ‘production. Because 
of the seasonal cy cles it was not possible to estimate which 
acid or which process gave more uniform results. 
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Moisture, 
Number 
of 





Laboratory 
Pilot Plant 


a 


Bulk 
density, 


g./ml. min. 


The Canadian Journal of Chemical Engineering, April, 1960 


2 day 
Gardner 
consistency, 


Stability, % 
Stir 
time, 


7 days 28 days 


12 36 
11 39 





The Upward Vertical Flow of 
Air-Water Mixtures 


III. Effect of Gas Phase Density on Flow Pattern, Holdup 


and Pressure Drop' 


R. A. S. BROWN?, G. A. SULLIVAN® and 
G. W. GOVIER* 


Results of measurements of pressure drop and 
holdup are reported over a range of air and water 
rates for the upward vertical flow of air-water mix- 
tures in a 1.50 inch I.D. tube with average air 
densities ranging from 0.092 to 0.552 lb./ft.? Super- 
ficial water velocities were between 0.0695 and 7.35 
ft./sec.. 

The data are analyzed by the method first sug- 
gested by Govier, Radford and Dunn and later used 
by Govier and Short. This involves the division of 
the flow range into four regimes on the basis of the 
pressure drop curves (to aid correlation of flow 
pattern and holdup data) and the separation of the 
unit pressure drop into hydrostatic and irreversible 
components. A superficial friction factor is calculated 
from the irreversible component of the unit pressure 
drop. 

The average density of the gas phase has a more 
or less marked effect on all three regime transitions. 
The transitions all shift similarly to lower air-water 
discharge volume ratios with increasing gas phase 
density. While the flow patterns were not directly 
observed in the investigation here reported, the flow 
pattern transitions would be expected to behave 
similarly. 

Over the range investigated, the gas phase density 
has little or no effect on the superficial friction factor. 
The holdup ratio is unaffected by the gas phase 
density in Regimes I and II but is significantly 
affected in Regimes III and IV. 


.o paper is the third of a series based upon an experimental 
study of upward vertical flow of gas-liquid mixtures. The 
object of the study is the development of a workable correlation 
for the effects on the two-phase flow of: 

(a) the flow rates of the two phases; 

(b) the diameter of the tube, and 

(c) the physical properties of the two fluids. 
The pertinent physical properties are believed to be density, vis- 
cosity and interfacial tension. The quantities which have been 
chosen as descriptive of the two-phase flow are: 

(a) unit pressure drop (i.e. pressure drop per unit length of 

tube) ; 
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(b) holdup ratio (i.e. the ratio of the volume ratio of the 
phases in the discharge mixture to that in the flow sec- 
tion), and 

(c) the flow pattern and its relationship to the unit pressure 
drop curve (i.e. the unit pressure drop versus the volume 
ratio of the phases in the discharge mixture). 

The volume ratios referred to above are in all cases the volume 
ratio of the gas phase to the liquid phase. 

The significance of the major variables, the thermodynamic 
relationships controlling the pressure drop and the relationship 
between the holdup ratio and the relative velocities of the phases 
are dealt with in detail in Part I of the series. 

The literature on two-phase gas-liquid flow was reviewed 
briefly in Part I and has been covered more thoroughly in two 
recent bibliographic articles». Most of the work on the subject 
has been devoted to the accumulation of data and the develop- 
ment of practical correlations with little, if any, theoretical basis. 
Notable among the exceptions to this generalization is the work 
of Calvert and Williams on annular flow. The work reported 
in this series was undertaken for the purpose of obtaining a com- 
prehensive set of self consistent data covering a wide range of 
the variables and to serve as a basis for the better physical 
understanding of the flow phenomena. 

The first two parts of this series have dealt, as does this, 
with the air-water system. Part I“) of the series was concerned 
with the effect of relative flow rates in a smooth bore tube hav- 
ing an inside diameter of 1.025 inches. For an average test 
section pressure of 36.0 p.s.i.a. and an average temperature of 
approximately 70°F. data were reported for discharge air-water 
ratios ranging from 0 to 248 and water rates varying from 
0.00040 to 0.0421 ft.*/sec. (superficial water velocities between 
0.0698 and 7.35 ft./sec.). Part II‘ presented the results obtain- 
ed in a study of the effect of diameter on the flow. Measure- 
ments were made under similar temperature, pressure and flow 
rate conditions with three additional tubes ranging from 0.63 to 
2.50 inches diameter. The work fteported here is primarily 
concerned with the effect of gas phase density on the two-phase 
flow, and is from a more general study of two-fluid systems. 
Investigations now in progress deal with the effects of viscosity 
and interfacial tension. 

The previously reported“: analysis of the flow system has 
shown that the overall unit pressure drop (here expressed in 
terms of feet of water per foot of height) may be separated into 
two components—that due to the hydrostatic head and that due 
to the irreversibilities associated with the flow. From the 
irreversible component of the pressure drop, a friction factor, 
fis may be defined based upon the Fanning equation and the 
superficial water velocity. “:© 
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:ental Equipment and Procedure 


xperimental equipment used to obtain the data reported 
here was essentially the same as that used in the previous 
investigs ‘ions of the system“:®, The test section was a 22.88 
foot length of 1.50 inch I.D. seamless copper tubing similar to 
that described by Govier and Short.‘ In order to obtain suffi- 
cient variation in air pressure to give a significant gas phase 
density change, an additional supply was obtained by parallel 
connection of two 125 SCFM air compressors with a maximum 
pressure of 200 p.s.i.g. This permitted operation at flowing 
pressures to 110 p.s.i.a. Tests at 18 and 36 p.s.i.a. were made 
using the laboratory air supply. 

The method of operation of the equipment was similar to 
that used previously.‘ Measurements were made at three gas 


Experi 
The . 


phase densities with the air-water system, corresponding to mid- 
point pressures of 18, 72 and 110 p.s.i.a. Air-water volume 
ratios ranged from 0.72 to 669 for water rates from 0.000851 to 
0.091 ft.*/sec. (superficial water velocities from 0.0695 to 7.35 
ft./sec.). In order to simplify the analysis of the data, the tests 
were made at approximately the same superficial water veloc- 
ities as were used previously.“:» As it was not practical to 
control the ambient temperature, some variation of flow temper- 
ature between tests was experienced, but it is not believed to be 
significant. 


Experimental Results 


The basic experimental data for the air-water system are 
tabulated in Table 1 (partial) for gas phase densities varying 


TABLE 1 (PARTIAL) 
EXPERIMENTAL DATA 


TuBE DIAMETER 0.125 FT. 
TEst SECTION LENGTH 22.88 FT. 
TEst SECTION Mip PoInT PRESSURE AS INDICATED 
TEstT SECTION Mip PoINT TEMPERATURE AS INDICATED 
SPECIFIC VOLUME OF LIQUID PHASE = 0.01603 FT.*/LB. 
SPECIFIC VOLUME OF SATURATED GAS PHASE AS INDICATED 
SUPERFICIAL WATER VELOCITIES APPROXIMATELY CONSTANT AT 0.868 FT./SEC. 








2 3 4 5 
Test 
Section 
Mid Point 


Pressure 
Psia 


Mid 
Point 
Temper- 
ature °F 


Saturated 
Air Rate 
ft.3/sec.* 


Water 
Rate 
ft.3/sec. 


Test 
Number 


Discharge 
Air-Water 
Ratio** 


6 


7 8 


Specific 

Volume 
Saturated 
Gas Phase* 

ft.3/lb. | 


Pressure 
Drop 
Regime 


| Pressure 
Drop 
ft. H.0/ft. 


Test 
Section 
Air-Water 


Holdup 
Ratio*** 





0.0107 0.0765 

0.174 
~ 155 
.022 
.0375 
.0519 


115 


27A 
28A 
29A 
30A 
31A 
32A 
33A 


345S 
346S 
350S 
351S 
352S 
353S 
354S 
355S 
356S 
357S 
358S 
359S 
360S 


121A 
122A 
123A 
124A 
125A 
126A 
127A 
128A 
129A 
130A 
131A 


186A 
187A 
188A 
189A 
190A 
191A 
192A 
193A 
194A 
195A 


18.0 


396 
214 
282 
170 
128 




















. . 
At test section temperature and average pressure. 


ia Column 4 divided by column 3. 
“*Column 6 divided by column 7. 
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from 0.092 to 0.552 Ib./ft.3 and for a water rate of 0.0107 
ft.3/sec. corresponding to a superficial water velocity of 0.868 
ft./sec. The complete tabular data for the same gas phase den- 
sities and for water rates ranging from 0.000851 to 0.091 ft.3/sec. 

are given in Table 1 1 deposited as Document No. 6163 with the 
American Documentation Institute Auxiliary Publication 
Project, Photo Duplication Service, Library of Congress, 
Washington 25, D.C. 

Referring to Table | (partial), column 2 gives the mid- -point 
pressure, column 3 gives the water rate in cu. ft./sec. and column 
4 gives the calculated flow rate of w ater-saturated air at mid- 
point conditions. The mid-point temperature is given in column 
5. The quotient of column 4 and column 3, the volume ratio of 
the air to water rates is given in column 6, while the volume 
ratio of air to water in the test section is given in column 7. 
The specific volume of the saturated gas phase at mid- -point 
conditions is given in column 8. The holdup ratio, (i.e. the ratio 
of column 6 to column 7), is given in column 9. The total 
pressure drop over the test section, expressed i in terms of feet of 
water per foot of height, is given in column 10. The pressure 
drop regime“: is shown in column 11. 


Interpretation and Correlation of Results 
Effect of Gas Phase Density on Flow Pattern 


It has previously been shown":® that the flow patterns 
occurring in the upward vertical flow of air-water mixtures 
are closely related to the pressure drop curve. In the earlier 
investigations, the flow was divided into four regimes, the bound- 
aries of which were taken, by definition, as the volume ratios 
corresponding with the two minima and the maximum appearing 
in the pressure drop curve. It was also shown that the flow 
patterns observed were related to these arbitrarily-defined 
regimes. 

The effect of a change in gas phase density is best illustrated 
by considering the pressure drop curves obtained with different 
gas phase densities at the same water velocity. Figure 1 shows 
the effect of gas phase density for a water velocity of 0.868 
ft./sec. Multiple ordinate scales are employed on this figure to 
avoid overlapping of the curves. The curves shown in the figure 
illustrate the trend of the transition points to lower volume 
ratios with increasing density of the gas phase. Because of the 
interrelation between the flow patterns and the flow regimes 
previously demonstrated“: this trend applies also to transitions 
from one flow pattern to another. 

Correlations of the loci of the regime and flow pattern transi- 
tions are shown in Figures 2 to 4. Figure 2 indicates that the 
locus of the first minimum of the pressure drop curve (the 
transition between Regimes 1 and Il) for all diameters and gas 
phase densities studied is represented fairly well by a single line 
on coordinates V, and R,p'/8D*’, The dependence on tube 
diameter shown here is the same as has been noted previously. ‘® 
The boundary lines between the different flow pattern regions 
are estimated from previous work. The locus of the second 
minimum of the pressure drop curve is correlated in Figure 4 
where V; is plotted against R,p'*D~**, The greater influence 
of tube diameter shown here as compared with previously pub- 
lished results‘® has arisen from a re-evaluation of the original 
data together with the results of additional experiments. 

As can be seen in the lists of symbols accompanying each of 
these figures, much of the previously published data was used 
along w vith the new data in developing the correlations. 


Effect of Gas Phase Density on Holdup 

The overall effect of the density of the gas phase on the 
holdup curves is shown in Figure 5 for a superficial water veloc- 
ity of 0.295 ft./sec. Multiple ordinates are employed on this 
figure also to avoid overlapping of the curves. ‘The holdup 
decreases rapidly with increasing gas phase density at higher 
values of R, and the maximum of the curve shifts to lower values 
of R,. At lower values of R,, in Regimes | and II, there is no 
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apparent effect of gas phase density. The holdup data have been 
correlated separately for each regime. 

Figure 6 shows the correlation of the holdup ratio for 
Regime I. The data form a broad band but show a segregation 
based upon superficial gas phase velocity. The upper portion of 
the band includes the bulk of the data at superficial gas phase 
velocities below 1.7 ft./sec. (lower limit of data is 1.3 ft./sec.). 
The data for gas phase velocities of 2.0 ft./sec. and higher lie in 
the lower part of the band. 


The Canadian Journal of Chemical Engineering, April, 1960 


density, ib/tt? 
0.092 


REGIME TI 


HOLDUP RATIO, He 
(see individual scales at right) 


y 


PRESSURE DROP 


+e a i 
UK 


50 
Ry +! 


Figure 5—Typical holdup ratio 


| 
| 
| 
+} 


| | 
REGIME IL 


average gos- 
phose density, | 


ibsft? 


0092 
0.184 
0 362 
0552 





symbol 








HOLDUP RATIO, Hp 


oii 





emeneeenetinsiaiiidiemusie 














x 
S 





REGIME IZ | 


HOLDUP RATIO, Hp 


1 
uv | 

overage gos phose 

density, ib/ft® fhe deme 


oos2 0164 0362 








SUPERFICIAL WATER VELOCITY, V_ 
Figure 9—Holdup ratio correlation—Regime IV. 


(ft/sec) 


In Regime II also the holdup ratio falls within fairly well 
defined limits when plotted against RyD'*. The data, which 
are shown in Figure 7, exhibit little systematic scattering. The 
holdup ratio is not affected significantly by the gas phase density 
in either Regime I or Regime II. 

The gas phase density does have an effect on the holdup ratio 
in Regimes III and IV however. Figure 8 shows the correlation 
for Regime III with HzD~° plotted against Ry with parameter 
Vp. The correlation of the holdup ratio in Regime IV is 
shown in Figure 9. These data indicate that, at constant super- 
ficial water velocity, Hx is inversely proportional to p°-®. 
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Effect of Gas Phase Density on Pressure Drop 

The friction factor is obtained from the observed pressure 
drop through Equations (1) and (2) previously reported.“ As 
was found for the data at a mid-point pressure of 36.0 p.s.i.a., 
fairly good correlation is obtained if the friction factor is plotted 
against DV, with a parameter of D*-°V; as shown in Figure 10. 
The abscissa is proportional to the superficial Reynolds Number 
of the water since the density and viscosity of the water are 
essentially the same for all tests. 

In the development of the correlation shown in Figure 10, 
all of the data for the air-water system were considered, includ- 
ing the results previously presented. “> The points shown on the 
figure are representative of the complete range of variables on 
which the correlation is based, and are not solely from the pres- 
ent investigation. The number appended to each point is the 
experimental value of the parameter, V¢D*-5, for the point. The 
data show that the density of the gas phase has no significant 
effect on the friction factor. The analysis of the additional data, 
however, has resulted in some minor adjustments in the chart 
previously presented.‘ It should also be noted that the depend- 
ence of the friction factor on tube diameter has now been confirm- 
ed for'a gas phase density of 0.092 lb./ft.’, the original determin- 
ation having been made for a gas phase density of 0.184 Ib./ft.’. 
The conventional single phase friction factor has been included 
in Figure 10 as the curve for a parameter value of zero. The 
dashed lines shown for low parameter values are approximate 
only and were obtained by interpolation. 

While the data used to develop this correlation’ covered the 
four flow regimes, it was not possible to get equally good 
representation in all of these. The correlation shown in Figure 
10 is applicable to flow in the first three regimes but gives low 
values of friction factor in Regime IV. The use of the correl- 
ation in Regime IV is not recommended. 


Conclusions 

The correlation presented in this paper should permit the 
prediction of flow pattern, holdup and pressure drop for the flow 
of immiscible gas and liquid water phases with average gas phase 
densities in the range 0.092 Ib./ft. — 0.552 Ib./ft. and tube 
diameters from 0.5 to 3.0inches. The correlations are restricted 
to cases in which the fractional change in gas density is small 
enough to justify the use of an average value. 





Nomenclature 
D = tube diameter, ft. 
gD { AF 
; Sg nas 
fi. == superficial friction factor = Oe 
2V,? \AX],. 
AF ol a é ; 
—]} = irreversibility effect, feet of water per foot of height. 
AX}, 
& = dimensional conversion factor, lb. ft./lb.psec.? 
G = flow rate of gas phase, Ib./nr. 
L = flow rate of water, lb./hr. 
AP = pressure drop, lb./ft.? 
R, = volume ratio of gas phase to water, Gu¢/Lv, 
Uv, = specific volume of water, ft.3/Ib. 
Uc = specific volume of gas phase, ft.3/Ib. 
V, = superficial water velocity based upon the total tube 
cross-section, ft./sec. 
Ve = superficial gas phase velocity based upon the total 
tube cross-section, ft./sec. 
AX = incremental height, ft. 
; 1 
p = density of gas phase = —, lb./ft.3 
Ug 
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